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INTRODUCTION 


This supplement contains information about the two main computer pro- 
grams (Langley computer program AOS 13, loading program, and Langley com- 
puter program A0457, mean camber program) used to obtain the results pre- 
sented in NASA Technical Note D-4427 along with two supplementary pro- 
grams (Langley computer program A1590, aspect ratio program, and Langley 
computer program A1591, pivot determining program) used in obtaining input 
data for them. 

In part I of this supplement, the input and some output variables for each 
of the two main programs (A0313 and A0457) are presented and pertinent com- 
ments made. Further, sample listings of input and output data are shown and 
the entire computer program listings provided. 

In part H of this supplement, a discussion of the two supplementary pro- 
grams (A1590 and A1591) is presented with a list of input data required, and 
sample input and output listings as well as program listings are given. 




Figure L- General planform for which lifting surface program is applicable. The x means that in the listing of geometry data the 
x- and y- locations of these points are found by dividing by the ratio of initial to final outer panel semispan. 



I. MAIN PROGRAMS (A0313 AND A0457) 


LOAD DISTRIBUTION PROGRAM A0313 


First card format: 6F12.5: 


(a) AR 


b2 


overall area 


Symbols for Input Data 


(overall area excludes the rhombus containing 


chord- extension) 

(b) CHI leading- edge inboard sweep, positive for sweepback, degrees 

(c) ALAMD leading-edge outboard sweep, positive for sweepback, 

degrees 

(d) PSI trailing- edge inboard sweep, positive for sweepback, 

degrees 

(e) BIRAT leading- edge break (distance from plane of symmetry to 

leading- edge break divided by semispan) 

(f) B2RAT trailing- edge break (distance from plane of symmetry to 

trailing-edge break divided by semispan) 

Second card format: 6F12.5: 


(g) TAPER 


(overall) 


tip chord 
root chord 


(root chord excludes chord- extension) j 


(h) DELT 


(i) XP 


(j) YP 

(k) CHDEXT 


(1) MACH 


change in the variable-sweep wing outer-panel sweep angle, 
positive if wing is swept rearward (zero for fixed wings), 
degrees 

x-location of the pivot point from the half root chord divided 
by the semispan when the tip is streamwise, includes the 
effect of the chord- extension 
y- location of the pivot point from the plane of symmetry 
divided by the semispan 

chord-extension (distance from rhombus leading-edge 
intersection with the plane of symmetry to the rhombus 
trailing edge divided by semispan) - if chord- extension 
is not 0, use only AJTEST = 2; however, if chord- 
extension is 0, values of AJTEST of 1 or 2 may be used 
(see fig. 1) ^ 

Mach number, must be less than 1 and should in practice be 
kept less than 0.9 


Tip stream- 
wise at 
DELT = 0 


Tip stream- 
wise at 
DELT = 0 


Tip stream- 
wise at 
DELT = 0 
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Third card format: 5F6.0, F6.2, F6.0: 


(m) CASE 

(n) SYM 

(o) CSTA 

(p) SSSTA 


case number 

symmetry code for loadings - if code is 1, the loadings are symmet- 
rical about the plane of symmetry; if code is 2, the loadings are 
antisymmetrical 

number of chordwise pressure modes and chordwise control points 
( a maximum number of 10 can be used) 
number of spanwise stations on a panel where the chordwise pressure 
modes are specified to act, includes the station at the plane of symmetry 
(a maximum number of 21 can be used); an equation which can be used 
as a guide in the determination of SSSTA is 


SSSTA 



Maximum combination value of CSTA and SSSTA are given in the following 
table: 


CSTA 

SSSTA 

1 

21 

2 

21 

3 

21 

4 

21 

5 

20 

6 

16 

7 

14 

8 

12 

9 

11 

10 

10 


(q) AJTEST if AJTEST = 1, the reference chord and reference area will be based on 

the total wing planform; if AJTEST = 2, the reference area and reference 
chord will be based on a wing planform which is determined by extending 
the leading and trailing edge of the outboard panel to the plane of sym- 
metry; if AJTEST = 2 and a variable-sweep wing is used with DELT 
other than 0, the reference wing will be obtained from that outboard panel 
when DELT is 0; note that AJTEST is not chosen indiscriminately, but 
should be chosen in connection with the value of the chord- extension as 
discussed earlier 
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(r) CLDESG desired lift coefficient at which the local loadings are to be calculated 

(s) TWADCM twist and camber code - set equal to 0 if the wing is flat and equal to 1 if 

the wing is warped; if set equal to 1, the local slopes of the control points 
when the root is at an angle of attack of 0 must be provided as input data; 
they are to be determined in the following manner: 

a l ~ ”(dx) + “twist 

due to camber 

and become the program terms CONST(JK,2); these terms, each of which 
is associated with a control point, are read in along each constant (x/c) 
row of control points starting nearest the leading edge and from root to 
right wing tip; the format for these numbers is 8F9.5 

MEAN CAMBER SURFACE PROGRAM A0457 

The geometric input data for this program is the same as for the program A0313 
except that the items AJTEST, CLDESG, and TWADCM are not used. The coefficients of 
the chordal loading function QP(J,N) are read in after the geometry items for each chord- 
wise pressure mode (starting with the cot ^ and ending with sin(CST A - 1)3) from the 
plane of symmetry to the right wing tip. The format for QP(J,N) is 8F9.5. 

COMMENTS ON INPUTS FOR PROGRAMS A0313 AND A0457 

The input data are based on the planform in a streamwise tip position with no inboard 
trailing- edge chord-extension. If the planform is of a fixed wing with unbroken leading 
and trailing edges and a skewed tip, in order to use these programs it is necessary to 
(a) determine its aspect ratio, (b) use the pivot program to find a pivot location for which 
in a lower sweep position the tip will be streamwise, and .(c) put the required geometric 
input data in with DELT set equal to the tip skew angle. This will result in the original 
fixed wing. 

A variable-sweep wing with trailing- edge chord- extension may with increasing 
sweep angle completely cover up the extension and intersect the fixed portion of the 
trailing edge, inboard of the chord- extension. If this is foreseen before that case is run, 
compute another wing for input whose aspect ratio and trailing- edge break are found by 
extending (spanwise) the fixed inboard trailing edge until it intersects the trailing edge 
of the outboard panel when the tip is streamwise. This new wing will therefore have 
no trailing-edge chord- extension, a different area, aspect ratio, and trailing-edge break 
value. 
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DEFINITIONS AND COMMENTS ON OUTPUT LISTINGS 


Some of the items of the output listings are not fully self-explanatory; therefore 
they are defined here. Also, some general comments are included where necessary. 

Geometry data : 

Perimeter points 
X LE REF 

Aerodynamic data (for 

CDl/CLA**2 
CDH/CLA**2 

LOCAL CIRCU 

SPAN LOAD ADD 

BASIC LOAD NO LIFT 


SPAN LOAD TOTAL 

The terms CENTER/PRESS-ADD, BASIC, and TOTAL refer to the chordwise loca- 
tion of the centroid of the various combinations of pressure loadings from the half root 
chord; the terms LOCAL/A. C.- ADD, BASIC, and TOTAL, which are obtained from the 
CENTER/PRESS terms, represent the distance from the centroid of the loadings to the 
local leading edge divided by the local chord and are listed in fractions of the local chord. 
Note that the local chords used are based on the "rounded" planform whose half-chords 
are given in the aerodynamic section of the listing of program answers. 

COMMENTS ON AERODYNAMIC DATA FOR PROGRAM A0313 

For variable- sweep wings, if the reference dimensions are taken to be those of the 
outer panel extended to the root, they are scaled internally to be in the same proportions 
to the wing in its new sweep position when divided by its new semispan as it was before 
the sweep change was made. These dimensions are used internally as scaled quantities, 
but are listed out in terms of the wing when the tip was streamwise. This means that the 


A0313): 


see figure 1 

location of the leading edge of the reference chord 
with respect to the overall half root chord 

due to additional loading only 
this has a sweep correction presented by Garner 
(ref. 13 of TN), also due to additional loading only 


= (W 

\ 2b ) a= \ 


rad 


T CLDESG 




, c L c a v/ 

0!=1 rad 


\ C 3lVL 


basic 


(— ) 
\ c av/ 


at ^L, desired 
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overall coefficients presented are based on the reference dimensions as listed, and 
the aerodynamic center is scaled appropriately so that it may be used directly with 
aerodynamic-center values at other sweep positions in forming plots of aerodynamic 
center as a function of sweep angle. However, if one desires to compute C ma from the 
chordwise location of the center of pressure, the reference chord value listed out and the 
moment reference point must be divided by the ratio of new-to-old semispan to bring the 
necessary quantities to the same scaling. Also, if one wants to compute the span load 
coefficient from the chord loading, the reference area must be divided by the square of 
the ratio of the two semispans when used in conjunction with the given lift coefficient as 
shown in the following equation: 

c l c _ 2 (chord load)/t>A=o\^ 

^L c av ^L®ref 

COMMENT FOR THE MEAN CAMBER SURFACE DATA OF PROGRAM A0457 

Note that the z/c terms in program A0457 contain the effects not only of camber 
but also twist and angle of attack. 

ADDITIONAL COMMENTS FOR POSSIBLE USE OF PROGRAM A0313 

Program A0313 can also be used to find the stability derivatives C^ and Cmq- 

To find C^p: 

(a) Set SYM = 2. and TWADCM = 0. 

(b) Replace DO 2 JK = 1, JKMAX 

2 CONST (JK,1) = 4.0 
With DO 2 JK = 1, JKMAX, 8 

2 READ (5,515) CONST(JK,l), CONST(JK+l,l), . . . CONST(JK+7,l) 

(c) Read in four times the linear twist in radians. This twist distribution is chosen 

to represent the variation of local angle of attack across the wing span which 
occurs due to a rolling velocity. A positive normal velocity is attained when 
the right wing tip is rolling up. 

(d) Take the number that appears at CROLL and divide by the rolling rate specified 

in terms of radians /second. The resulting number is 

To find Cmq^ 

(a) Set SYM =1.0. and TWADCM = 0. 
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(b) Same as in (b) for C i . 

Jr 

(c) Read in four times the linear camber in radians. This camber distribution is 

chosen to represent the variation in local angle of attack across the wing 
chord which occurs due to a pitching velocity. A positive normal velocity is 
attained when the nose is pitching up and a zero normal velocity occurs at 
either the center of gravity or c/4, whichever is selected for the wing to pitch j 
about. j 

(d) Take the number that appears at CMA and divide by the pitching rate specified 

in terms of radians/second and multiply by (2/CREF). The resulting number 

is (-Cm q ). 
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LOAD DISTRIBUTION PROGRAM A0313 




-p 

m 

U 

•H 

P4 


H 

Qi 


m 


sa 

O ft •<*> 
o g w 
0) S c3 

CQ CQ O 


rt 

IS 

Q 

! 

HH 

<D 


Cli 

03 


-d- 



CO to CM 

cm oo en 


m m m 

O 



N OK\ 

hcd 

tH 

co r-^ o 

in 



tOHJ 

O HlPi 

rH 

-d 00 CM 

m o 



o m 

00 oo 

-d 

tn o tH 

o to 



HHO 

O iH o 

O 

o o o 

• • 


• • 

• • • 

0 0# 

0 

0 0# 

o o 


o o 

1 1 1 

III 

8 

4- 1 




en N H 

d O J- 

CO 

■d* H rH 




N Oto 

m cm 

CM 

in -d OJ 




oens 

tn o o 

OO 

tH m rH 

en 


m 

o m cm 

oo 

m 

tO »H H 

r^. 


m 

i — 1 tH O 

o HO 

O 

o o o 

eg 


CO 

* » • 

• * * 

0 

0 0 0 

en 


o 

1 I 1 

l i + 

1 

4- 1 

CM o 


CM 





• • 


.# 0 





o o 


o 

to o o 

en NN 

tH 

into s 




LnNOO 

IDOJ- 

en 

m -d oo 




-=Mn o 

o ento 

-d 

oo co en 




en m cm 

s cno 

m 

tO rH O 




O tH O 

o oo 

O 

O o o 

tn 



I 3 1 

1 1 + 

i 

Hh 8 








in 







00 



m to 

o in o 

CO 

en in o 

o o 



rH oo en 

n ho 

o 

rH 00 O 

• B 


• » 

oo *h m 

in nm 

«H 

in m co 

CO O 

€> 

o o « 

co tn h 

to en h 

m 

N CM O 

in 

o 

tH 

O rH O 

o o o 

o 

o o o 




1 V 1 

1 8 4* 

8 

+ 1 

en 


CM 

CM CM CO 

03 OH 

m 

en r*^ tH 

rH O 

CM 

tn to m 

tn in 

co 

in in to 

O 

• 

• 

HNOO 

en -d-cn 

to 

en in 

O 

rH 


CO CM -d- 

tn CD rH 

CM 

•d\CM O 

O LA 



O iH »H 

o oo 

o 

o o o 

* • 


• * 

• # » 

• 00 

0 

0 0 0 

tn o 

« 

CM O • 

I 9 I 

1 1 4* 

8 

8 1 

»H 

rH 

CO tH 

tH tO 

en i^ cm 

in 

r- to eg 




co en in 

in co cm 

iH 

N LO N 




cm r*"* 

m rH to 

eg 

to tn cm 


» 

• 

cm 

in CD CM 

CM 

-d m o 



O 

OHH 

o oo 

o 

o o o 



CM 

0 0* 
1 1 K 

0 0 0 
t 1 4- 

l 

0 0 0 
8 8 

O O 






e « 


• * 





o tn 

* 

m o * 

o en cm 

to om 

o 

m o in 

co 

CO .d* 

O CO CM 

nino 

tH 

•d e- to 




co to 

r^v co m 

h- 

in rH O 




to rH -d* 

-d- com 

rH 

-d O 




OHH 

o o o 

O 

O O o 

en in 

« 

O-O- * 

8 8 g’ 

8 6 + 

1 

1 

•d* CM 

tH 

CO N H 





cc m 


O CO 





to 


CD Jf 

to o tn 

m 

rH 

to CM oo 

tH tH 


-d* O 

h tn m 

en tooo 


in eg -d- 

* * 


• * 

rH CM Ut 

O -d“ 0 

tH 

noa d 

in o 


rH * 

lOrl d 

-d* com 

Hddo 


in 

O 

O «H H 

o oo 

o 

o o o 



O 

»»0 

0 « 0 

0 

0 0 0 



CM 

8 1 8 

8 8 + 

8 

I 



10 




11 





13 



Program Listing 
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noon 


PROGRAM L IF SURF t INPUT , OUTPUT, TAPE5= INPUT, TAPE6=0UTPUT ) 
MODIFIED LIFTING SURFACE THEORY 

FOR VARIABLE SWEEP PLANFORMS 


DIMENSION Y(4l) , ETA (41) ,C(41),D(4l) ,8 (4 1 ) , IP IVO ( 100 ) ,SUML ( 100, 100 ) 
1 » PH I ( 10 ) » VE ( 41 ) ,CONST( 100,2) ,DF( 10) ,AL( 10,41) ,SA( 41 ,41) ,Q( 10,41 ) ,Q 
2T (10,41 ) , CHDLD(41)» CHDLDT (41 ) »CHDLTQ( 41 ) ,CLL ( 41 ) ,CL B( 41 ) ,CLLTOT ( 41 
3 ) , GAMMA ( 41 ) ,GAMMAT( 41 ) ,G AMT MAC 41 ) ,AMU(4 1 ) ,AMUT( 41 > , AMUTOT (4 1 ) , Z I CP 
4(41) ,ZICTP(41> ,ZICTPT(41) ,SLC(41) ,BLDT( 41 ) , BSPLDT (41) ,TAW( 10) , JJJ( 
510) , INDEX (100,2 ) , ALPHA! (41) 

REAL MACH ,LOAC( 4 l ) , LOACTO( 4 1 ) »LOACT (41) 

EXTERNAL F0FT1 , F0FT2, F0FT3, F0FT4,FUFT5,F0FT6 ,F0FT7, F0FT8,F0FT9, FOF 
IT10 

COMMON XSUB ,YSUB ,ETASUB ,CSUB, DSU3 

500 FORM AT ( 6F12 .5) 

501 FORMAT ( 1H1, 58X, 13HGE0METRY DATA///) 

502 FORMAT ( 1H1, 57X, 16HAER0DYNAM IC DATA///) 

503 FORMAT ( 50X, 5HX2 PP=, F9.5 ,5X, 5HY2PP=, F9 .5 ) 

504 F0RMAT(50X» 5HX3NP=, F9 .5 ,5X, 5HY3NP=» F9.5 ) 

505 FORMAT ( 50X, 5HX4 PP=» F9.5 ,5X, 5HY4PP— » F9»5 ) 

506 FORMAT ( 50 X, 5HX5 PP=, F9.5 » 5X, 5HY5PP='» F9.5 ) 

507 FORMATt 50X, 5HX6NP=, F9 .5 ,5X, 5HY6NP=, F9.5 ) 

508 F0RMAT(50X,5HX7PP=,F9.5,5X,5HY7PP=,F9.5) 

509 FORMAT ( 40X12HCASE NUMBER=,F6.0, 5X, 14HSYMMETRY CODE= ,F5.0,5X, 12HMAC 
1H NUMBER=,F9.5) 

510 FORMATt 1H0) 

511 FORMATt IX »I 2 , 10F9.5 , 10X , 2F9 • 5 ) 

512 FORMAT ( 14X, 3HADD,6X , 5HBAS IC , 4X5HT0T AL ,4X, 3HADD, 6X , 5HBASIC,4X, 5HT 
10TAL ,4X ,3HADD, 6X,5HBASIC ,4X , 5HT0TAL, 5X , 3HADD ,5 X, 7HN0 LIFT,2X,5HTO 
2TAL ) 

513 FORMAT ( 23X, 7HN0 LIFT,1X,9HAT CLOESG , 10X , 7HN0 LIFT,1X,9HAT CLDESG,1 
10X, 7HN0 LIFT , IX »9HAT CLDESG ,18X,9HAT CLDESG) 

514 FORMATt 15X,97HIF SYMMETRY CODE IS EQUAL TO l, THE SPAN LOADING IS S 
1YMMETRICAL, OTHER THAN 1,IT IS ANT IS YMMETR ICAL ) 

515 FORMATt 8F9. 5 ) 

516 FORMAT ( 55X, 42HCH0RDAL LOAD FACTORS, Q, FOR ADDITIONAL LOAD) 

517 FORMAT { 1X4H SPAN2X,4H2Y/B3X, 5HCH0RD, 4X ,5HCH0RD,4X,5HCH0RD, 4X ,6HCENT 
1ER,3X,6HCENTER,3X,6HCENTER,3X,5HL0CAL,4X,5HL0CAL,4X,5HL0CAL,5X,4HS 
2PAN ,4X» 5HBASIC»4X »4HSPAN) 

518 FORMAT ( IX ,4HSTA. »9X , 4HL0AD, 5X ,4HL0AD, 5X ,4HL0AD, 5X» 5HPRESS ,4 X,5HPRE 
1SS, 4X,5HPRESS ,4X,4HA»C. , 5 X, 4HA.C . , 5X, 4HA.C. ,6X, 4HL0AD,4X,4HL0AD, 5X 
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2 »4HL0AD ) 

519 FORMAT! 26X, 4HSPAN ) 

520 FORMAT! 1 X ,4HSPAN2X» 4H2Y/B3X » 5HLQCAL »4X, 5HL0CAL, 4X.5HL0C AL,4X,5HL0C 
1 AL , 4X, 5HL0C AL »4X.» 5H LOCAL » 4X , 5HL0C AL ,4X, 5HL0CAL , 4X , 5HL0C AL * 14X, 5HL0 
2C AL 4X , 5HL0C AL ) 

521 FORMAT i l X ,4 HSTA . , 9X ,5HC IRCU , 4X, 5HC I RCU, 4X, 5HC IRCU,4X, 4HLI FT ,5X ,4H.L 
1IFT,5X,4HLI FT,5X,5HPITCH,4X,5HPITCH,4X»5HPITCH, 14X, 4HHALF , 5X,6HL0C 
2ATN) 

522 FORMAT! 30X13HASPECT RAT I0=, F9.5 , 5X, 14HPLANF0RM AREA*, F9. 5 , 5 X, 14HAV 
1ERAGE CH0RD=,F9.5// 10X21HMEAN GEOMETRIC CHQRD= , F9 . 5 , 5X, 59HX LOCATI 
2 ON OF THE LEADING EDGE OF THE MEAN GEOMETRIC CH0RD*»F9.5/ 30X39HY 
3L0C AT ION OF THE MEAN GEOMETRIC CHORD*, F9, 5//) ' 

523 FORMAT! ///50X,60HCH0RDAL LOAD FACTORS, QT» FOR THE LOAD DUE TO TWIST 
1 AND CAMBER) 

524 FORMAT! 15X, 11HR00T CHORD* ,F 9. 5 , 5X , 1 OHTI P CHORD* ,F9. 5*5X, 15HF0REWIN 
1G C HORD* »F9.5»5X»20 HOVER ALL TAPER RAT 10= , F9 . 5/ / 45X1 1HY LE BREAK*, F 
29.5 » 5X» 1 1HY TE BREAK* »F9. 5//25X17HX PIVOT LOCAT ION* ,F9. 5 , 5X 1THY PI 
3V0T LOC AT ION* * F9 . 5, 5X » l 7HZ PIVOT LOCATI ON* , F9. 5 //55X, 19HTE CHORD E 
4 XT E NS ION* »F 9. 5 ) 

525 FORMAT! 5F6. 0,F6. 2, F6.0) 

526 FORMAT! ////40X59HL0CATI0N OF PERIMETER POINTS FOR THE PLANFORM USE 
ID AS INPUT) 

527 FORMAT I//5X35HNUMBER OF CHORDWISE PRESSURE MODE S= ,F 5.0 , 5X72HNUMBER 

1 OF STATIONS SPANWISE ON A PANEL WHERE PRESSURE MODES ARE DEFINED* 

2 * F5 • 0) 

528 FORMAT! 5X ,4 1HCH0RDW I SE LOCATION OF CENTER OF PRESSURE*, JF9, 5 , 5X, 60H 
1A.C. IN FRACTION OF CREF MEASURED FROM LEADING EDGE OF CREF=,F9.5/ 
26X , 40HSPANW I SE LOCATION OF CENTER OF PRESSURE* , F9 . 5 ) 

529 FORMAT! 50X,9HX1PP= 0. 1 OX ,9HY1PP= 0.) 

530 FORMAT! 57X, 12HCASE NUMBER*, F6.0) 

531 FORMAT! 10X, L5HCLA PER DEGRE E* , F9 . 5 , 5X , 1 5HCMA PER DEGREE* , F9 .5 » 5X, 1 
1 1 HC 0 1 / C L A ** 2= , F 9 . 5 , 5X, 12HCD 1 1 /CL A** 2*»F9. 5/// ) 

532 FORMAT ( 52 X, 7HXI* 0.,10X,7HY1= 0.) 

533 FORMAT! IX , I 2 , 4F9 . 5, F8. 5 , El 1 . 3 ,F8. 5 , F7 .5 * El 1 . 3 ,4F9. 5 ) 

534 FORMAT { 26X, 14 , 2X , 10F9 .5 ) 

535 FORMAT !52X,3HX2=»F9.5»5X,3HY2=»F9.5 ) 

536 F0RMAT(52X,3HX3*,F9.5,5X,3HY3=,F9.5) 

537 FORMAT! 52X, 3HX4* ,F9.5,5X,3HY4*,F9.5 ) 

538 FORMAT! 14X , 3HADD ,6X , 5H8 AS IC ,4X, 5HT0TAL, 4X , 5HC0EFF ,4X, 5HC0EFF ,4X » 5H 
l COEFF , 3X ,8H COEF F, LE , IX, 8HC0EFF, LE , IX, 8HC0EFF ,LE , 12X , 5HCH0RD ,4X, 6HM 
2 I DC HD ) 

539 FORMAT! 25X5H ST A . ,5X , 14 , 5X , I 4, 5X , 14 ,5X, I 4, 5X, 14 ,5X , 14, 5X, I 4, 5X, 14, 
1 5X, 14, 5X, 1 4 ) 
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540 FORMAT ! 50X, 5HX6 AP=» F9.5 ,5X, 5HY6AP=, F9.5 ) 

541 FORMAT ( 5X7HREF AR=, F9.5 ,2X, 14HREF AREA, SR£F-»F9 .5,2 X, 17H HOME NT REF 
1 PQINT=,F9.5,2X15HREF CHORD, CREF= ,F9.5» 2X,9HX LE R£F=,F9.5) 

542 FORMAT! 31X, 1HN, 6X,2HLS,9X,4H2Y/810X ,3HX/C,7X,8HDELTA CP , 5X, 8HDELTA 
1 CP5X, 8H0EL TA C P/76 X, 3HADD, 9X, 5 HBAS IC , 8 X, 5HT0TAL/87X , 7HN0 LIFT»5X, 
29HAT CLDESG) 

543 FORMAT ( 28X, 14,3 X , 14 ,6X,F9. 5 , 5X, F9.5 ,4X, F9 . 5 , 4X, F9. 5 ,4X, F9 .5 ) 

544 FORMAT C 2E16 .8 ) 

545 FORMAT (52X»3HX5=»F9.5,5X,3HY5=,F9.5 > 

546 FORMAT ( 2X ,4 HCL A= , F9 . 5 , 2X »4HCMA= , F9. 5, 2X , 6HCR0LL =,F9 .5, 3X , 20HCL , TWI 
1ST AND CAMB ER=»F9.5,2X»20HCM» TWIST AND CAMBER=, F9.5,2X, 10HCL,D£SIG 
2N=» F9.5 ) 

547 FORMAT ( 23X, 7HN0 LIFT,1X,9HAT CLDESG , IX, 3HADD , 6X , 5HB ASIC , 4X, 5HT0TAL 
1 , 4X , 3HA DD,6X,5HBASIC,4X , 5HT0TAL/50X , 7HN0 L£FT,1X,9HAT CLDESG, 10X, 7 
2HN0 LIFT , IX »9HAT CLDESG) 

548 FORMAT ( 5X ,2 3HLE INBOARD SWEEP ANGLE = , F9 .5»5X,32HLE INITIAL OUTBOAR 
ID SWEEP ANGLE—, F9»5,5X, 30HLE FINAL OUTBOARD SWEEP ANGLE= , F9 .5/5X , 
223HTE INBOARD SWEEP ANGLE=, F9. 5 , 5X, 32HT E INITIAL OUTBOARD SWEEP AN 
3GLE-,F9.5 ,5X,30HTE FINAL OUTBOARD SWEEP ANGLE=, F9 .5//5X,40 
4HCHANGE IN OUTER PANEL SWEEP ANGLE , DELT A= ,F9. 5 , 4X , 26HP I VOT CANT AN 
5GLE IN PITCH=,F3.0,4X,25HPI VOT CANT ANGLE IN ROLL= , F3 .0/7 ) 

549 FORMAT ( 5 IX, 4HX6 A= , F9 . 5, 4X, 4HY6A= ,F9 . 5 ) 

550 FORMAT (52X, 3HX6= , F9 . 5 , 5X, 3HY6= , F9.5 ) 

551 FORMAT l 2X , L 22HTDT AL WING PLANFORMi MEAN GEOMETRIC CHORD AND ITS LOC 
l ATI ON ARE INVALID IF THE PLANFORM HAS A TRAILING EDGE CHORD EXTENS 
2I0N) /// ) 

552 FORMAT { //// 55X , 20HREFERENCE DIMENSIONS//) 

553 FORMAT ( ///3 3X ,72HL0CAT ION OF PERIMETER POINTS FOR PLANFORM TO BE U 
1SED IN THE COMPUT AT IQNS/40X , 57HWHEN NONDIMENSIONALI ZED BY THE SEMI 
2 SPAN RATIO GIVEN ABOVE) 

554 FORMAT { //50 X, 26HFQR THE ADDITIONAL LOADING) 

555 FORMAT! 40X,57HWHERE THE ORIGIN IS AT THE HALF ROOT CHORD, PQSI T I VE 
IX AFT//) 

556 FORMAT! 52 X, 3HX7= , F9 . 5 , 5X , 3HY7= , F9.5 ) 

557 FORMAT ( 4E 16 . 8 ) 

558 FORMAT! 25X, 33HI SEMI SPAN AT FINAL OUTBOARD SWEEP ,37H/SEMISPAN AT IN 
1ITIAL OUTBOARD SWEEP) =,F9.5 ) 

559 FORMAT! 10X, 31HANGLE FOR ZERO LIFT, ALPHA ZER0=,F9. 5, 15X, 44HP ITCHING 
1 MOMENT COEFF. AT ZERO LIFT, CM Z£RO= ,F9 . 5/30X , 53HR00T BENDING MOME 
2NT COEFFICIENT AT ZERO L I FT ,CBMROOT = , F9 .5 ) 

INPUT DATA 
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c 

C IF THE WING DOES NOT HAVE A LEADING EDGE BREAK SET BIRAT AND CHI 

C EQUAL TO ZERO 

C 
C 

C IF THE WING DOES NOT HAVE TRAILING EDGE BRE'AK SET B2RAT AND PSI 

C EQUAL TO ZERO 

C 
C 

C IF THE WING DOES NOT HAVE A VARIABLE SWEEP OUTER PANEL SET OELT , 

C XP AND YP EQUAL TO ZERO 

C 

C 

C IF THE SPAN LOADING IS TO BE SYMMETRICAL SET THE SYM CODE EQUAL TQ 

C ONE, IF ANT I SYMMETRICAL SET THE SYM CODE EQUAL TO TWO 

C 
C 

C IF THE WING HAS A TRAILING CHORD EXTENSION SET AJTEST EQUAL TO 

C TWO. THIS WILL BASE ALL OF THE AERODYNAMIC COEFFICIENTS ON THE 

C PLANE OR M OF THE OUTER PANEL EXTENDED TO THE ROOT 

C 
C 

C SET CLDESG EQUAL TO ONE UNLESS A PARTICULAR VALUE OF THE LIFT 

C COEFFICIENT WOULD BE USEFUL IN SCALING THE LOADING DISTRIBUTION 

C 
C 

C IF THE WING HAS TWIST AND/OR CAMBER SET THE TWADCM CODE EQUAL TO 

C ONE. OTHERWISE SET IT EQUAL TO ZERO 

C 
C 
C 

1 READ! 5, 500) AR» CHI , ALAMO , PS I , BIRAT, B2RAT 
I FI EOF, 5) 10,5 

5 READ I 5, 500) TAP ER,DELT, XP ,YP»CHDEXT , MACH 

READl 5, 525) CASE ,SYM,C STA, SSST A , AJ TEST , CLDESG, TWADCM 


VARIBLE SWEEP GEOMETRY PROGRAM 


PI=3. 14159265 
RAD-180. /PI 
Q8AR= 1.00 
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B0= 1.0 
YMI N=l. 0 
YMAX=1.0 
YMAX1=1.0 
SIGM=0. 

RHH=0. 

ZP=0. 

CHDEXX= CHDE XT 

B£TA=SQRT(1.-MACH**2> 

I SYM=SYM 

JMAX=CST A 

I SS ST=SSSTA 

NMA X=2* I SSS T— 1 

NNI 1= I SSST+ 1 

JKMAX=JMAX* ISSST 

JTE ST=AJTEST 

NTWACM=TWA0CM+1. 


CHI=CHI /RAD 
AL AMD=ALAMD/RAD 
PSI=PSI /RAD 
DELTA=DELT/RAD 
SIGMA=SIGM/RAD 
RH0=RHH/RlD 
TANC=TAN(CHI ) 

TANL=TAN( AL AMD) 

TANP=T AN ( PS I ) 

TANDE=T AN(P I /2. +DELTA ) 

B1=B1RAT*B0 

82=B2RAT*B0 

CR=B0*(4./AR-82RAT*TANP-TANC*(B1RAT*(B1RAT-B2RAT-1. ) )-TANL* 

1131 RAT* ( B2R AT-BlRAT+1 . ) -B2R AT ))*(!•/( B2RAT* ( 1.— TAPER ) + { l. + TAPER)) ) 
0MEGA=ATAN < ( 1. / ( 1. -B2R AT ) ) *{ (TAPER-1.) *(CR/B0) +B1RAT*( TANC-TANL) 
1— 82RAT*TANP+TANL) ) 

TAN0=TAN (OMEGA) 

EOMEG=OMEGA+DELtA 
ALAME=AL AMD+DELTA 
TAME=T AN ( EOMEG) 

TANA=T AN ( AL AME) 

CR=CR+CHDEX T 
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X2=-CR/2. 

Y2=0.00 

Z2=0.00 

X3=-CR/2.+Bl*TANC 

Y3=B l 

Z3=0.00 

X4=-CR/ 2. +3 1* { TANC- TANL )+8Q*TANL 

Y4=B0 

Z4= 0. 00 

X5= CR/2 . +B2* l TAN P— T ANO ) + B0*T AND— C HD EXT 

Y5=B0 

7.5=0.00 

X6A=CR/2.+B2*TANP 

Y6A=B2 

Z6A=0 .00 

X6= X6A-CHDE XT 

Y 6= Y6 A 

Z6=0. 00 

X7=CR/2 . 

Y7=0. 00 
Z 7=0. 00 
X2PP=X2 
Y2P P=Y2 
Z2PP=Z2 
X6A P=X6 A 
Y6AP=Y6 A 
Z6AP=Z6 A 
X 7PP=X7 
Y7PP=Y7 
Z7PP=Z7 

IF( DELTA. EQ.O.) GO TO 360 

X1PP=XP+(X3-XP) *COS ( S IGMA )*COS (CELT A) +{ Y3-YP ) *CQS 1 RHO ) *SIN I DELTA) 
L + ( Z 3— ZP ) * ( COS ( S IGMA )*SI Nt RHO ) *S INI DELTA )— SIN l SIGMA) *COS { RHO) ) 
Y3PP=YP — (X3-XP) *COS l SIGMA )*SIN( CELT A) +( Y3-YP ) *COS( RHO) *COS( DELTA) 
1+ { 13- IP ) * (COS ( S IGMA ) * SIN(RHO)*COS( OELT A) +5 IN (S IGMA ) *SIN< DELTA ) ) 
Z3PP=(X3-XP )*SIN (SIGMA) *COS ( RHO)-( Y3-YP )*SIN< RHO) 

1 + ( Z3-ZP ) *C0 S { S I GMA ) *COS iRHO ) + ZP 

X4PP=XP+(X4-XP) *COS IS IGMA) *COS( DELTA) +( Y4-YP)*C0S I RHO) *.$ IN (DELTA) 
1 + IZ4-ZP )*(COS IS IGMA) *S I N(RHO)*S INI DELTA) -SINT SIGMA )*COS( RHO)) 

Y4P P= YP- ( X4 -XP) *COS (SIGMA)*SIN { DELT A) +( Y4-YP ) *C OS ( RHO ) *COS { DELTA ) 
I+(Z4-ZP)*IC0S< SIGMA)* SI NCRHO )*CGSl DELTA) *S INI S IGMA )*SIN1 DELTA) ) 
Z4P P=(X4— XP )*SIN(SIGMA) *COS (RH0)-(Y4-YP )*SIN(RHO) 

1+1 Z4-ZP ) *C0 S( SI GMA) *COS ( RHO )+ZP 
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X5P P=XP+ ( X5— XP) *COS (SIGMA ) *COS( DELT A) + ( Y5-YP ) *CQS( RHO )*S IN ( DELTA) 
1+(Z5— ZP)*(COS(S IGMA )*SI N(RHO )*SI NIDELTA )— SINI S IGMA ) *COS( RHO ) ) 
Y5PP=YP— (X5-XP) *COS ( S IGMA )*SI N ( DELT A) +1 Y5-YP)*COS( RHO )*COS( DELTA) 
H-IZ5-ZP)* (COS (SIGMA)* S IN(RHO ) *COS( DELT A) +S IN (S IGMA )*SIN( DELTA ) ) 
Z5PP=(X5— XP ) *SI N{ S IGMA ) *COS {RHO)-CY5-YP )*SIN{RH0) 

1 + ( Z5-ZP ) *CQ S< S I GMA ) *COS ( RHO ) +ZP 

X6PP=XP+(X6-XP)*C0S( SIGMA )*COS( DELTA )+( Y6-YP )*COS( RHO) *S INI DELTA) 
1 + (Z 6— ZP ) * (COS ( SIGMA ) *SI N( RHO ) *S I NIDELTA )— SIN (SIGMA) *COS l RHO ) ) 

Y6P P=YP— ( X6— XP) *COS I S IGMA )*SIN( DELTA) +( Y6-YP )*C OS (RHO) *COSC DELTA) 
1+IZ6-ZP)* (COSTS IGMA)* Si N (RHO) *COS( DELTA )+SIN( S IGMA )*SIN( DELTA ) ) 
Z6P P= ( X6- XP ) * SI N ( S I GMA ) *COS ( RHO ) - ( Y6-YP ) *S I N ( RHO ) 
l + ( Z6-ZP ) *CO S ( S I GMA ) *COS ( RHO ) +ZP 
IF ( ( Y3-Y2PP ) . EQ .0. ) GO TO 208 
A2PP=(X3 -X2PP ) /( Y3 -Y2PP) 

GO TO 209 

208 A2PP=0. 

209 IF ( ( Y6— Y7PP ) .EQ.O. ) GO TO 210 
A6PP=(X6 — X7PP) / (Y6 -Y7PP) 

GO TO 211 

210 A6PP=0. 

211 IF ( (Y4PP-Y3PP).EQ.0. ) GO TO 212 
A3PP=( X4PP-X3PP ) / (Y4PP-Y3PP ) 

GO TO 213 

212 A3PP=0. 

213 IF( ( Y5PP— Y6PP) .EQ.O* ) GO TO 214 
A5P P= ( X5P.P-X6PP ) / ( Y5PP-Y6PP ) 

GO TO 215 

214 A5PP=0. 

215 IF ( ((X3-X2PP)**2+(Y3-Y2PP)**2).EQ.O.) GO TO 216 
G2PP= ( Z3 -Z2PP ) / ( ( X3 — X2PP ) **2+( Y3 - Y2PP ) **2 )**. 5 
GO TO 217 

216 G2PP=0. 

217 IF( l ( X7PP-X6 ) **2+( Y7PP-Y6 )**2 ) • EQ.O» ) GO TO 218 

G6PP= (Z6 -Z7PP ) / ( ( X7PP— X6 )**2+(Y7PP-Y6 )**2)**.5 

GO TO 219 

218 G6PP=0. 

219 IF{ (A3PP— A2PP). EQ.O. ) GO TO 220 

X3NP= ( A3PP*X2PP-A2PP*X4PP+A2PP*A3PP*(Y4PP-Y2PP) ) / ( A3PP-A2PP) 
Y3NP= ( X2PP-X4P P+A3PP*Y4PP-A2PP*Y2PP) / ( A3PP-A2PP) 

GO TO 221 

220 X3NP=X2PP 
Y3NP=Y2PP 

221 Z3NP-( (X3NP-X2PP)**2+(Y3NP-Y2PP)**2 )**. 5*G2PP+Z2PP 
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IF ( ( A6P P—A5 PP ) • EQ*0 * ) GO TO 222 

X6NP= ( A6PP*X5PP-A5PP*X7PP+A5PP*A6PP*iY7PP-Y5PP ))/l A6PP-A5PP) 
Y6NP={ X5PP- X7PP+A6P P*Y7P P— A5PP*Y5PP ) / { A6PP— A5PP) 

GO TO 223 

222 X6NP— X7 PP 
Y6N P-Y7 PP 

223 Z6NP=( ( X6NP-X7PP}**2+(Y6NP~Y7PP)**2)***5*G6PP+Z7PP 
IFCCHDEXT.EQ.O. ) GO TO 342 

X6NP=X6PP+{ Y6-Y6PP) *A5PP 
Y6NP=Y6 

Z6NP= ( X6AP- X6NP ) * ( Z6A-Z6 ) / ( X6A-X6 ) +Z6AP 

342 IFIY4PP.GE* Y5PP) GO TO 411 
IFC Y4PP.LT. Y5PP) YMAX=Y5PP 
YMAX1=YMAX 
GO TO 410 

411 YMAX=Y4PP 
YMAX1-YMAX 

410 X2P P=X2 PP/YMAX 
Y2P P-Y2 PP/Y MAX 
Z2PP=Z2PP/YMAX 
X3N P-X3NP/YMAX 
Y3N P-Y3NP/Y MA X 
Z3NP=Z3NP/Y MAX 
X4PP=X4PP/YMAX 
Z4P P=Z4 PP/YMAX 
X5PP-X5PP/Y MAX 
Z5PP= Z5 PP/YMAX 
X6AP = X6 AP/Y MAX 
Y6AP = Y6 AP /Y MAX 
Z6AP=Z6 AP/Y MAX 
X6NP=X6 NP/Y MAX 
Y6NP= Y6 NP/YMAX 
Z6NP= Z6 NP/Y MAX 
X7P P = X7 PP/YMAX 
Y7PP=Y7 PP/YMAX 
Z7PP=Z7PP/YMAX 
XP= \P/ Y MAX 
YP= YP/ Y MAX 
ZP= ZP/YMAX 
CR-CR/Y MAX 
Y4P P = Y4PP/YMAX 
Y5P P= Y5 PP/Y MAX 
CHO£XT=CHDE XT/Y MAX 



IF < ABSI ALAME-CHI KG£. 0*0001 74) GO TO 765 

X3NP=X2PP 

Y3NP=Y2PP 

Z3NP= Z2 PP 

765 IF(ABSCEOMEG-PSI). GE.O. 000174. OR.CHOEXT . NE. 0- 0) GO TO 767 
X6NP=X7 PP 

Y6NP=Y7PP 

Z6NP=Z7PP 

X6AP=X7PP 

Y6AP=Y7PP 

Z6AP-Z7PP 

767 IF ( CHDEXT .NE»0*0) GO TO 766 
X6AP-X7PP 
Y6AP-Y7PP 
Z6AP=Z7 PP 

766 IFf Y6AP.EQ* Y6NP) GO TO 1050 
CH0EXX=0.0 

GO TO L 051 

1050 CHDEXX— X6AP— X6NP 

1051 SAOO=CHDEXX *Y6NP- 
YMA X=l* 00 

YHI N-AM INli Y4PP , Y5P P ) 

B1RAP=Y3NP 

B2RAP=Y6NP 

IF i Y6NP-Y7PP) 362,361,361 
362 X6NP1=X6NP 
Y6NPi~Y6NP 
Z6NP1-Z6NP 

X6NP-=X6NP1“Y6NP 1*TANE 
Y6NP-=0 : * 0000 000 

TANO={ Z6NPI-Z5PP I / C ( X6NPI-X5PP ) **2+ 1 Y6NP1-Y5PP ) **2) *#* 5 
Z6NP- Z6NPI-C IX6NP— X6NP 1 ) **2+( Y6NP*^Y6NP 1 ) **'2 )** *3*1 ANU 
CR= X6NP—X2P P 
OR I GNN- f X6N P-X7PP ) / 2* 

X2PP^X2PP*-0RIGNN 
X3NP=X3NP~0RIGNN 
X4P P*.X4 P P-0 R I GNN 
X5P P==X5 PP—O R I GNN 
X6NP=X6NP-QRIGNN 
X6A P = X6 A P—0 R I GN N 
X7P P = 0* 00 
32RAP=0.00 
TANP=0* 000000 
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GO TO 361 

360 X3NP-X3 
Y3NP=Y3 
Z3NP=Z3 
X4PP=X4 
Y4PP=Y4 
Z4PP=Z4 
X5PP=X5 
Y5PP=Y5 
Z5P P=Z5 
X6NP=X6 
Y6NP=Y6 
Z6NP=Z6 

SADD=CHDEXX*Y6NP 

81RAP=Y3NP 

B2R AP=Y6NP 

361 TANOP=TANE-TANP 
TANOL=T ANE-TANA 
TANOC=T ANE-T ANC 
TANLC=TANA-TANC 
TANLP^T ANA-TANP 
TANPC=TANP-TANC 
IFtY4PP.GE.Y5PP) XCBLEM= 

1 (CR+Y3NP*TANLC-Y6NP*TAN0P+Y5PP*TAN0U*t Y4PP*l-CR/2.-Y3NP*TANLC) 

2 + (Y4PP+Y5PP )/2.*l Y4PP*TANA-l-CR/2.-Y3NP*TANLC) )-(Y4PP**2 
3+Y4PP*Y 5PP+ Y5PP *#2 ) # TANA/ 3. ) 

IF ( Y4PP .LT.Y5PP) XCBLEM= 

1 ( CR +Y3NP*TANLC— Y6NP*T ANOP+Y 4PP*T ANOL ) *( Y5PP*(-CR/2.-Y3NP*TANLC 
2+Y4PP* ( T ANA-TANDE) )+(Y5PP+Y4PP)/2.*(Y5PP*T ANDE- {-CR/2.-Y3NP*TANLG 
3+Y4PP*( TANA-TANDE) ) )-l./3.* ( Y5PP**2+Y5P P*Y4PP+Y4PP**2 ) *TANOE ) 
CR=CR— CHOEXX 

S=2 .* (- Y3NP**2*TANLC/ 2. +Y6NP**2*T AN0P/2 . +YMIN* ( CR+Y 3NP*TANLC 
l-Y6NP*TAN0P+YMIN*TAN0L/2. ) + SADO + 

2 t YMAX-YMI N) /2 .* ( CR+Y3NP*T ANLC-Y6NP*TAN0P+YM IN*T ANOL ) ) 

CR=CR+CHDEXX 


C8AR « XC BAR , AND YCBAR ARE NOT VALID FOR A DISCONTINUOUS TRAILING EDGE 


YCBAR= 2 ./S*( CR/2.* YMIN ** 2 +Y 6 NP **3*T ANaP/6.-Y3NP**3*TANLC/6. 
1+YMIN**2*< YMIN/3.*TAN0L+Y3NP/2.*TAN LC-Y6NP/2.*TAN0P ) 

2 + 1 YMAX+2.*YMIN)*( YMAX-YMI N)/6.*(CR+Y3NP*TANLC-Y6NP*TAN0P 
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3+YM I N*T ANOL ) 1 
CAV=S/(2.*YMAX> 

IF{ B1RAP-B2RAP) 303,304,304 

303 CBAR=2» /S*( CR**2*Y3NP-»-CR*Y3NP**2*TANPC+ Y3NP**3*TANPC**2/3. + 

l { CR-*-Y3NP*TANLC ) *#2* ( Y6NP— Y3NP )— ( Y6NP**2— Y3NP**2 ) * ICR+Y 3NP*T ANLC ) * 
2TANLP+iY6NP**3-Y3NP**3)*TANLP**2/3. + 

3 C CR + Y3NP*TANLC— Y6NP*T ANOP )**2*l YM IN— Y6NP ) + ( YMI N**2-Y6NP**2 ) *TANOL* 
4 (CR+ Y3NP*TANLC— Y6NP*TAN0P ) + ( YMI N#*3— Y 6N P**3 ) *T ANOL **2/3 • + 
5(YMAX-YMIN)/3.*( CR+Y3NP*TANLC-Y6NP*TAN0P+YMIN*T ANOL )**2 ) 

XCB AR=2 . /S* (-YMI N/2 . *CR**2+CR/2. * ( 1 . 5*Y3NP**2*T ANLC— Y6NP**2*TANQP 
l/2.+YMIN**2*(3.*TAN A— TANE ) / 2.-3 . 0*Y MI N* Y3NP*TANLC+YMI N* Y6NP* 
2TAN0P) + Y3NP **3*T ANLC*I2.*T AN A— 4« *TANC+TANP) /6. + Y6NP**3*TANA*T ANOP / 
36.+YMIN**3*TANA*TAN0L/3.-Y3NP**2*YMIN*TANLC**2+(YMIN*Y3NP*Y6NP 

4- Y3NP/2.*Y6NP**2)*TANOP*TANLC+Y3NP/2.*YMIN**2*TANLC*l2.*TANA-TANE) 

5- Y6NP/2 . *YM I N**2*TANA*TAN0P+ 

6XCBLEM ) 

GO TO 301 

304 CBAR=2. /S*( CR**2*Y6NP+CR*Y6NP**2*TANPC+ Y6NP**3*TANPC**2/3. +(CR 

1— Y6 NP*T ANOP ) **2* { Y3NP— Y6NP ) + t CR— Y6NP*TAN0P ) * ( Y 3NP**2— Y6NP**2) 
2*TAN0C+ ( Y3NP**3-Y6NP**3 )*TAN0C**2/3. + 

3(CR+Y3NP*TANLC-Y6NP*TAN0P)**2*( YMIN-Y3NP) +(CR+Y3NP*TANLC-Y6NP* 
4TAN0P ) *TANQL*( YMIN**2— Y3NP**2) + l YMIN** 3-Y3NP**3) *TANQL**2/3. + 

5 ( YMAX-YMIN ) /3 .* ( CR+ Y3NP*T ANLC-Y6NP*TANQP+YMI N*TANOL )**2 ) 
XCSAR=2./S*(-Y6NP/2.*CR**2+Y6NP**2/2.*I CR*TANC-CR/2.*TANPC)+Y6NP** 
13/3.*TANC*TANPC- { Y3NP-Y6NP )*CR/2 . *( CR-Y6 NP*T ANOP )+ { Y 3NP**2- 

2Y6NP**2 )/2. *(-CR*TAN0C/2.+TANC*(CR-Y6NP*TAN0P) )+{Y3NP**3-Y6NP**3) / 
33.*(TANC*TANOC) + < YM I N-Y3NP) * ( ICR+Y3NP*TANLC-Y6NP*TAN0P) *(-CR/2 

4.-Y3NP*TANLC) ) + ( YMI N**2-Y3N P**2 ) / 2. * ( TANA* l CR+Y 3NP*TANLC-Y6NP*TAN0 
5P)+TANOL*<-CR/2.-Y3NP*TANLC ) )+( YMIN**3- Y3NP**3) /3. *TANA*TANOL+ 
6XCBLEM ) 

301 ARN=4.*YMAX**2/S 
ARB=ARN*B£T A 
CR=CR/BET A 


THE FLAT PLATE ANGLE OF ATT ACK , CONST ( JK , 1 ) , I S SET EQAUL TO ONE 
RADIAN 


DO 2 JK=1 , J KMAX 
2 CONSTUK, 11=4.00000000 
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THE LOCAL ANGLE OF ATTACK OUE TO TWIST AND/OR CAMBER, CONST! JK,2» , 
ARE EITHER READ IN OR ASSIGNED A ZERO VALUE HERE OEPENOING ON THE 
VALUE OF THE TWADCM CODE 


IF! NTWACM.EQ.l) GO TO 1002 

READ (5, 515) I CONST! JK,2 ) , JK=1, JKMAX ) 

DO 8 1=1, JKMAX 
8 CONSTII * 2 )=4. *CONST ( 1,2) 

GO TO 1004 

1002 DO 1003 JK=1, JKMAX 

1003 CONST! JK,2)=0.0 
C 

1004 CHI 8= AT AN (T ANC/BETA ) 

ALAMB=ATAN! TANA/BETA) 

PS I 8=AT AN (T ANP/BETA ) 

OMEGB=AT AN! TANE/BETA) 

T ANCB=T AN ( CH I B) 

TANL8=T AN! ALAMB ) 

T AN P B=T AN ! P S I B ) 

TANOB=T AN(OMEGB) 

TANLCB=T ANL B-TANCB 
TANLPB=T ANL B-TANPB 
TANPCB=TANP B-TANCB 
TAN0CB=TAN08-TANC8 
TANOLB=T ANOB-TANLB 
TANOPB=TANO B-TANPB 
ETA! ISSST)=0.Q 
Y ( I SSST ) =0. 0 
DO 7 NP=NNI I , NMAX 
ANP=NP 

ETA(NP)=SIN! ( ANP-SSSTA)*PI /(2.0*SSSTA) ) 

Y(NP)=ETA(NP) 

CH0SUB=0. 

IF (ETA! NP) . GT.( Y6NP-0.10) . AND.ETA INP ) .L E. Y6NP. AND.CHDE XX.NE .0. ) 
1CHDSUB=CHDEXX*( l.-( Y6NP-ETA (NP) )/.10) /BETA 
IF (ETA(NP).GT.Y6NP. AND.CHDE XX.NE. 0. ) CHDSUB=CHDEXX/BETA 
IF! Y3NP.GE.Y6NP) GO TO 307 
IF(ETA(NP).GE.Y3NP) GO TO 309 
C(NP) = (CR+ETA(NP)*T ANPC B-CHDSUB) /2.0 
D( NP)=ET A (NP )*( TANP B+TANCB ) /2.0-CH0 SUB/2.0 
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DIFF3=ETA(NP)-0. 

312 IFl ITTU.NE.2) GO TO 314 
IUSX=NP-1 

-I UST=NP— 2 

C 1 I USX ) = 1 10 . *C { I USX )+2. *C I IUSD1/12. 

D 1 1 U S X ) = 1 1 0 . * 0 1 I USX ) + 2 . *D C I US T ) ) / 1 2 , 

ITTU=1 

314 IFIDIFF3.LT .O..OR.D IFF3 .GT. *01) GO TO 323 
I TTU=2 

323 IFl NP-NNI I ) 7,6,7 

309 IFl ETAl NP).GE.Y6NP) GO TO 327 

C 1 N P ) = ( CR +Y 3 NP * T ANL C B-E T A I N P ) *TANL P B-CH D S UB ) / 2 . 0 
D (NP )=(— Y3NP*TANLCB+ETA INP ) *1 TANPB+TANLB.)— CHDSUB)/2.0 
DIFF3=ETA(NP)— Y3NP 
GO TO 312 

327 IF 1 ETAl NP) . GE. YMI N) GO TO 331 

Cl NP )-( CR+Y3NP*TANLCB— Y6NP*T AN0PB+ETA1NP ) *TANOL 8-CHOSUB) /2. 0 
D(NP)=(-Y3NP*TANLCB-Y6NP*TAN0PB-»-ErAlNP)*(TAN08+TANLB)-CHDSUB 1/2.0 
DIFF3=E TAl NP)— Y6NP 
GO TO 312 

331 IFl YMIN.EQ.Y5PP) GO TO 340 

CATY4=CR+Y3NP*T ANLCB— Y6NP*T AN0PB + Y4PP*T ANOLB-CHDSUB 
DAT Y4=X4PP/BETA+CAT Y4/2. 

Cl NP ) = . 5*1 CATY4— CAT Y4*l ETAl NP )-YMIN )/ I YMAX-YMIN ) ) 
D(NP)=DATY4+(X5PP/BETA-DATY4) *1 ETAl NP )-YMIN)/( YMAX-YMIN) 
DIFF3=ETA(NP )-YMIN 
GO TO 312 

340 C(NP)=1 Y4PP-ETA (NP) ) / (2. 0*1 Y4PP-Y5PP) ) * (CR+Y3NP*TANLCB-Y6NP 
1*TAN0PB+Y5PP*T ANOLB-CHDSUB) 

DATY5=(-Y3NP*T ANLCB— Y6NP*TAN0PB+Y5P P*(TANLB + TANOB )— CHDSUB ) / 2. 0 
Dl NP )=DATY5 +1 X4PP/BET A-DATY5 ) *1 ETAINP )- YMIN) / 1 YM'AX- YMIN) 
DIFF3=ETA(NP)-YMIN 
GO TO 312 

307 IFl ETAl NP) . GE.Y6NP) GO TO 311 

C 1 N P ) = 1 CR +E TA 1 N P ) *T ANPC 6-CHDSUB )/2.0 
D(NP)=ETA(NP)*( TANPB+TANCB) /2.0— CHDSUB/2.0 
DIFF6=ETA(NP)-0. 

313 IFl ITTU.NE.2) GO TO 316 
I USX=NP— 1 

IUST=NP— 2 

Cl IUSX)=(10.*C( I USX ) + 2.*C 1 1 UST ) ) / 12. 

DIIUSX)=(10.*D(I USX ) +2 . *0 1 1 UST 11/ 12 . 

I TTU=1 



316 IF( DIFF6.LT . 0..0R.DIFF6 .GT. .01) GO TO 325 
I TTU=2 

325 IF(NP-NNIl) 7,6,7 

311 I F ( E T A ( NP ) • GE . Y 3NP ) GO TO 329 

C ( N P ) - ( CR-Y 6 N P* T ANO P B+E T A { N P ) *T ANQC B-*CH DSUB } / 2 . 0 
D( NP ) =( ~Y6NP*TAN0PB+ETA ( NP ) *( TANOB+T ANCB) — CHDSUB) /2 .0 
DIF F6=ET AI NP )— Y6NP 
GO TO 313 

329 IF(ETAINP).GE.YMIN) GO TO 333 

C ( NP ) = { CR+Y 3NP*T ANLCB-Y6NP*TAN0P B+ETA (NP )■*( T AN0L8 )-CHDSUB 1/2. 0 
Q(NP)=(-Y3NP*TANLCB-Y6NP*TAN0PB+ETA{NP) *1 TAN0B+TANL8)-CHDSUB)/2.0 
DIFF6=ETAINP)— Y3NP 
GO TO 313 

333 IF C YMIN.EQ. Y5PP ) GO TO 346 

CATY4=CR+Y3NP*T ANLCB— Y6NP*TAN0PB+Y4PP*T ANOLB—CHOSUB 
DATY4=X4PP/BETA+CATY4/2. 

C (NP ) = . 5*(CATY4-CATY4*( ETA(NP)— YMIN) / (YMAX— YM'IN.J.) 
D(NP)=DATY4+(X5PP/BETA-DATY4)*(ETA( NP J-YMIN) / ( YMAX-YM IN J 
DIF F6=E TA ( N P) -Y M I N 
GO TO 313 

346 C(NP)=( Y4PP-ETA(NP» ) / ( 2.0*< Y4PP-Y5 PP! I *<CR+Y3NP*TANLC8-Y6NP 
1*TANQPB+Y5PP*TAN0LB— CHDSUB) 

DATY5=(— Y3NP*TANLCB— Y6NP*TANQPB+Y5PP*(T ANLB+TANOB)— CHDSUB )/ 2.0 
D ( NP ) =DATY5 + ( X4P P/BET A-DATY 5 ) *(ETA( NP J-Y.MIN)/ C YMAX-YM IN ) 
DIFF6=ETA(NP)-YMIN 
GO TO 313 

6 C( ISSST)=(5.*CR+2.*C(NNII ) ) /I 2. 

D( I SSST ) =0( NN 1 1 )/6. 

7 CONTINUE 

DO 187 N=NN 1 1 »NMAX 
KSU=NMAX+ 1— N 
ETA ( KSU )=— ETA( N ) 

Y( KSU ) =— Y ( N ) 

C(KSU)=C(N> 

D(K SU) = Dt N) 

187 CONTINUE 

CH I I=RAO*CH I 
ALAM=R AD* ALAMO 
ALA DL=RAD*A LAME 
PS I BB=R AD*P S I 
OME G=R A D*OM EGA 
EOMG=RAQ*EOMEG 
CR=CR*B ETA 
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CTB=(X5-X4) /YMAXI 

TAPER=CTB/CR 

CFB=Y3NP*TANC 


WRITE INPUT OAT A 


WRITE (6,501) 

WRI TE ( 6 , 509 ) CA SE , S YM, MACH 
WRITEI6 ,514) 

WRITE (6,527) CSTA,SSSTA 
WRI TE (6 ,510 ) 

WRITE (6, 526) 

WRI TE (6 ,555 ) 

WRI TE (6 ,532 ) 

WRIT£(6»535) X2,Y2 
WR I TE (6,536) X3,Y3 
WRI TE (6 ,537 ) X4,Y4 
WRI TE (6 , 545 ) X5,Y5 
WRI TE ( 6 ,550 ) X6,Y6 
WRI TE ( 6 » 549 ) X6A.Y6A 
WRITE (6, 556) X7,Y7 
WRI TE (6 , 510 ) 

WRI TE (6 ,558 ) YMAXI 
WRITE(6,510) 

WRI TE (6 ,553 ) 

WRI TE 16 ,555 ) 

WRI TE ( 6 , 529 ) 

WRI TE ( 6 ,503 ) X2PP.Y2PP 
WRI TE (6 ,504) X3NP,Y3NP 
WRI TE ( 6 ,505 ) X4PP»Y4PP 
WRITE(6,506) X5PP,Y5PP 
WRI TE(6 ,507 ) X6NP, Y6NP 
WRITE (6, 540) X6AP,Y6AP 
WRI TE ( 6 , 508 ) X7PP,Y7PP 
WRITE(6,510) 

WRI TE ( 6 , 551 ) 

WRI TE (6 ,522 ) ARN, S, CAV.CBAR , XCBAR ,YCBAR 

WRITE (6, 548) CH 1 1 , ALAM, AL AOL, PS I BB, OMEG , EOMG, OELT , S IGM, RHH 
WRITE (6, 524) CR , CTB ,CFB , TAPER , BIRAP , B2RAP ,XP ,YP , ZP.CHDEXX 
C 
C 
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TO DETERMINE WHAT REFERENCE DIMENSIONS TO USE 


IF ( JTES T. NE • 1 ) GO TO 400 

CREF=CB AR 

YR£F=YM AX 

SREF=S 

XLREF=XCBAR 

REFPT=XC8AR+CBAR/4. 

ARNN=ARN 
GO TO 401 

400 CRSW = CTB-(T ANO— TANL ) /YMAX1 
TAPRSW=CTB/CRSW 

CR£F= ( 2 . #CR SW/3 . ) *( l.+TAPRSW+TAPRSW**2) / ( 1. +TAPRSW) 

YREF=1. /YMAX1 
SREF=<CTB+CRSW) /YMAX1 

XLREF=2./SREF*(-CRSW**2/ (2.*YMAX1 )— GRSW / ( 2. *YMA XI **2 ) *1 ITANO-TANL 
l )/2.-TANL)+TANL/(3.*YMAXl**3)*{TAN0-TANL) ) + ( X4— Y4*TANL ) /YMAX1 
2+CRSW/2. 

REFPT=XLREF+CREF/4. 

ARNN= ( 2 . /YM AX 1 ) *#2/ SREF 


MAIN PROGRAM 


401 DO 3 NTR=l,41 
VE ( NTR) =0.0 

B (NTR) =0.0 

CHOLD ( NTR) =0 .0 
GAMMA (NTR)=0.0 
CLL (NTR) =0.0 
CHDLDT ( NTR) =0.0 
ZICP (NTR) =0.0 
LOAC ( NTR) =0.0 
SLC (NTR) =0.0 
BLOT (NTR) =0.0 
BSPLDT ( NTR) =0.0 
CLB ( NTR ) =0.0 
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CHDLTO(NTR)=0.0 
CLLTOT ( NTR)=0,0 
GAMMAT { NTR) =0.0 
GAMTMAf NTR) =0.0 
ZICTP ( NTR) =0.0 
LOACT (NTR) =0.0 
ZICTPT(NTR)=0.0 
LOACTO( NTR) =0.0 
AMU (NTR) =0.0 
AMUT (NTR) =0.0 
AMUTOT ( NTR) =0.0 
ALPHAI ( NTR) =0.0 
3 CONTINUE 

DO 590 NIR=l,4l 
DO 590 NTR=L,4I 
590 SA( NI R» NTR) =0.0 
DO 78 KS = l » JMAX 
PKS=KS 

PHI l KS)=( 2. 0*PKS*PI ) / (2. 0*CSTA+l.O) 

00 78 NU=I S SST» NMAX 

JK=(KS-2)*ISSST+NU+1 

ANU=NU 

DO 14 N= 1 » NMAX 
AN=N 

VE( N)=COS( ( (AN— SSSTA)*PI)/(2.0*SSSTA) ) 

SA( N,N)=4./ (2.*SSSTA)*VE(N) 

NNUD=I ABS (N— NU) 

IF ( NNUD.NE.O) GO TO 9 

B(N)=(2.0*SSSTA)/(4.0*C0S(( ( ANU-SSSTA)*PI )/ ( 2.0*SSSTA> ) ) 

GO TO 14 

9 IF( MOD ( NNUD >2).EQ.O) GO TO 12 

8(N)=VE (N)/ ( (2.0*SSSTA)*(ETA(N)-Y(NU) )**2) 

GO TO 14 
12 B(N ) =0. 0 
14 CONTINUE 

DO 589 NP=I SSST »NMAX 
NPNUD=I A8S( NP-NU) 

IF(NPNUD.EQ.O) GO TO 589 

IF l MOD ( NPNUO* 2 ) .EQ.O) GO TO 589 

SA(NU»NP)=2.0*B(NP)/SSSTA*C0S( ( ANU- $S ST A) *P I /( 2 .0*S SSTA ) ) 

IT=NMAX+1— NU 

ITT-NMAX+I-NP 

SAl NU » I TT ) = 2 . 0*6( IT T) /SSSTA*COS ( ( ANU-SS STA) *P I / ( 2.0*SSSTA) ) 



noon 


SAI IT »NP) =SAI NU 1 1 TT ) 

SA( IT,ITT)=SA(NU,NP) 

SAt NP,NU)=SA(NUfNP) 

SA( ITT, IT)=SA(NU,NP) 

589 CONTINUE 

DO 78 J = 1 * J MA X 
AJ=J 

DO 30 N=1,NMAX 

AK=0.0 

AN=N 

IFIN.NE.NU) GO TO 16 
IF(J.EQ.l) GO TO 18 
I F C J— 2 ) 20,19,20 

18 AK=2.0*PHI (KS1+2.0* SINIPHKKSI) 

GO TO 21 

19 AK= PH I ( KS . 5* SIN(2.0*PHI(KSn 
GO TO 21 

20 GA= (SINC <AJ-2.0)*( PHKKSI1 1 l/IAJ-2.0) 

AK=GA— C S I N ( ( AJ 1* (PHI (KS))))/AJ 

21 PARTL=8 (N) *AK 
A=0 . 0 

DO 25 NUP=1 , NMAX 
NUPNU=I ABS( NUP-NU) 

IFI NUPNU.EQ.O) GO TO 25 

I FI MODI NUPNU, 2) .EQ. 0 1 GO TO 25 

SSND=ABS(Y( NU)-ETA(NUP) 1 

IFI SSND.EQ.O.) GO TO 25 

ANUP=NUP 

UURR=ANUP— SSSTA 

A=A+ I ( COS ( ( UURR#P 11/(2. 0* SSSTA) 11**2) *ALOG( SSND) 

25 CONTINUE 

IFIJ.NE.il GO TO 28 

0F< 1)=— 1. 0/12,0*1 SIN MPH I (KS) 1/2.0) )*(.SINM PHI IKS) 1/2, 0)1) 

GO TO 29 

28 DF 1 J 1 = 1 AJ— 1 .0 1 * (COS ( (AJ — 1.0 1* ( PHI (K S 1 1) 1 

29 VL =1.0 /( ( ( C( N) 1**2 1 *2 «Q*SSSTA*VE( N) *S IN (PHI IKS) ) ) 

ALt J,N)=VL*DFI J )*( ( .25*SSSTA)*{ 1.0-2.0*( VE( N1 ** 21-ALOG ( 4. 0) l-Al+PA 
1RTL 

GO TO 30 


CHORDAL INTEGRATION SUBROUTINE 

SOLVES FOR THE CHORDAL INFLUENCE FUNCTION VALUES 
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16 XSU8=— C { NU) *COS IPHI ( K S ) )+D!NU) 

YSU8=Y { NU) 

ETASUB=ETAI N) 

CSUB— C( Ni 
DSUB=D!N) 

GO TO (351,352,353,354,356,357,358,359, 1070,1071) 
1071 CALL GAUSS! 0. ,PI , 3, SUM1 0 , F0FT10 ) 

AK=SUM10 
GO TO 355 

1070 CALL GAUSS! 0. ,PI ,3, SUM9,F0FT9) 

AK=SUM9 
GO TO 355 

359 CALL GAUSS! 0.,PI,3,SUM8 »FOFT 8 ) 

AK= SUM8 
GO TO 355 

358 CALL GAUSS! 0.,PI,3,SUM7,F0FT7) 

AK= SUM7 
GO TO 355 

357 CALL GAUSS! 0. ,PI ,2, SUM6 ,F0FT6) 

AK=SUM6 
GO TO 355 

356 CALL GAUSS! 0. ,PI ,2, SUM5 ,F0FT5) 

AK= SUM5 
GO TO 355 

354 CALL GAUSS! 0. , P I , 2, SUM4 ,F0FT4 ) 

AK=SUM4 
GO TO 355 

353 CALL GAUSS! 0.»PI,2» SUM3 , F0FT3 ) 

AK= SUM3 
GO TO 355 

352 CALL GAUSS! 0. »PI»2»SUM2»F0FT2) 

AK= SUM2 
GO TO 355 

351 CALL GAUSS! 0, ,PI,2, SUM1,F0FT1) 

AK=SUM1 


355 AL! J, N) =-B! N) *AK 
30 CONTINUE 

DO 78 NP=ISSST, NMAX 
I=! J-2) *1 SSST+NP+1 
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IF( NP.EQ.ISSST) GO TO 73 
NR=NMAX+1-NP 
IF ( ISYM.NE. I) GO TO 77 
SUML( JK, I ) = AL ( J * NP ) + AL ( J » NR ) 

GO TO 78 

77 SUMLI JK,I)=AL(J,NP)-AL( J,NR> 

GO TO 78 

73 IF ( ISYM.NE. I) GO TO 75 
SUML ( JK » I ) = AL ( J » NP) 

GO TO 78 

75 SUML(JK,I)=0. OOOOOOOO 

78 CONTINUE 

MATRIX SOLUTION SUBPROGRAM 

SOLVES FOR CHORDAL LOAD MODIFICATION FACTORS Q AND QT 


CALL MAT INV ( SUML,JKMAX , CONST, 2, DETERM, IPIVO, INDEX, 100, ISCALE) 


DO 821 J=I, JMAX 
DO 821 NP=I SSST , NMAX 
IXX=(J-2)*I SSST+NP+1 
QT( J ,NP )=CONST ( IXX,2)*QBAR 
821 Q( J ,NP) =CONST( IXX, i)*QBAR 


DU 760 N=1 , NMAX 
Cl N )=C( N) *B ETA 
760 D(N)=0(N)*BETA 
DO 439 J-l , JMAX 
00 150 NP= I SSST, NMAX 
IF( JMAX.GE.2) GO TO 751 
QI 2 »NP) =0.0 
QT ( 2 ,NP )=0. 0 

751 CHDLDI NP)=PI*IQI 1 »NP) +. 5*Q( 2,NP ) ) 

GAMMA(NP)= CHDLDI NP ) / (4 .0*QBAR*YMAX I 
CLL I NP ) =2.0*YMAX*GAMMAI NP)/C(NP) 
CHDLDTI NP ) = P I *( QT 1 1 , NP ) + . 5*QT 1 2 , NP ) ) 
150 CONTINUE 

DO 439 NP=NNI I , NMAX 

NR=NMAX+1— NP 

IF! ISYM.NE. 1) GO TO 152 
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Q! J »NR) =Q{ J » NP) 
QT! J,NR)=QT!J,NP> 
GO TO 439 

152 Q ( J > NR ) =— Q ! J » NP ) 
QT(J»NR)=QT { J,NP) 
439 CONTINUE 


DO 650 NP=NNII,NMAX 
NR= NM AX + 1 -N P 
IF! ISYM.NE. 1) GO TO 649 
GAMMA! NR) =GAMMA(NP) 
CHDLD! NR )=CHDLD ! NP) 
CHDLDT! NR)=CHDLDT !NP) 

GO TO 650 

649 GAMMAINR)=-GAMMA!NP ) 
CHDLD! NR )=-CHDLDt NP ) 
CHDLDT! NR) = CHDLDT! NP) 

650 CONTINUE 


BL=0.0 

AMX=0.0 

AMY=0.0 

BLT=0.0 

AMYT=0. 0 

DO 154 N=l» NMAX 

IE! JMAX.GE.2) GO TO 754 

Q12 ?N)=0.Q 

QT( 2»N) =0.0 

754 RU=PI*VE ( N) /! 2. *SSSTA) 

BL=BL+RU*CHDLD( N) 

AMX= AMX +RU* ET A ! N ) *C HOLD I N ) 

BLT =BLT+RU* CHDLDT (N) 

IF! JMAX.GE. 3) GO TO 877 

Q13»N)=0.0 

QT!3»N)=0.0 

877 AMY=AMY-RU*PI*! ! D(N )-XCBAR-CBAR/4« ) *( Q( 1 ,N)+.5*Q< 2, N) J— Cl N) *!.5*Q{ 
il,N)+.25*Q!3,N) ) ) 

AMY T=AMYT— RU*PI *! !D!N)-XCBAR-CBAR/4. ) *( QT ! 1 ,N )+*5*QT! 2,N) )-C(N)*!. 
15*QT!1*N)+.25*QT!3,N) )) 

154 CONTINUE 

DO 155 NP-=I SSST * NMAX 
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IF( ISYM.EQ.2.AND.NP.EQ.ISSST) ZICPU SSST)=D(NP ) 

IF( ISYM.EQ.2 . AND. NP» EQ. I S SST 1 GO TO 11 

ZICP(NP)=(-C(NP)/2.0)*( (Q(ltNP)+.5*Q(3,NP) ) /( Q 1 1 ,NP > + . 5*Q< 2 ,NP) ) ) 
l+DINP) 

11 LOAC(NP)* (ZICP(NP)— (Q4NP I— C (NPI ) ) /{ 2. G*C( NP ) ) 

SLCfNP) = PI*(Qil,NP)+.5*Q(2,NP))/(BL/I2.*YMAX)) 

IF { NTWACM.EQ. 1 ) GO TO 4 

BID T ( NP ) =CH Dl DT { NP ) / ( BLT/ ( 2 ,*YMAX ) ) 

GO TO 155 
4 BLOTC NP 1=0. 

155 CONTINUE 

AMYP=AMY-BL* I (XCBAR+CBAR/4. l-REFPT) 

AMYPT=AMYT-8LT* 1 1 XCBAR+CBAR/4. )— REFPT ) 

CM=AMYP/ { Q3AR*SREF*CREF ) 

CMT=CM/RAD 
CLT=BL/ (QBAR*SREF) 

CLTT=CLT/RAO 

CLTWST=BLT/(QBAR*SREF) 

CMTAC=AMYPT/(QBAR*SREF*CREF ) 

ALPHZO=— CLTWST/CLTT 
CMZERO=CMTAC+CMT*ALPHZO 
CRL=AMX/ ( QB AR#SREF*2 * *YREF ) 

ZICPT=— AMYP/BL+REFPT 
TOAC=(Z ICPT— XLREFI/CREF 


RATO=CLTWST/CLT 
RTTG=RATO-CLDESG/CLT 
DO 1000 NP= I SSST t NMAX 

BSP LOT ( NP )= i BLDT t NP ) — SLC ( NP ) ) *CLTWS T*SR EF/S 
BLDT ( NP ) =BSPLDT ( NP) +SLC ( NPT*CLD£SG*SREF/S 
CLB ( NP) =BSPLDT ( NP )*CAV/ ( 2.*Cl NP) ) 

IF( NTWACM.EQ. 1) GO TO 200 
CHDLDT(NP)=CLB(NP)*2.*C(NP)*QBAR 

ZICTP(NP)=(-C(NP)/2.Q)*t((QT(l»NP )+ .5*QT< 3, NP > l-TQC 1,NP )+. 5*QI 3 ,NP 
1) )*RATO)/(I QT(1,NP) +.5*QTI2,NP) )-(Q(l,NP)+.5*QI 2,NP) )*RATO) )+D(NP) 
ZICTPT<NP)=«-C<NP)/2.)*(((QT<1,NP)+.5*QT{3,NP))-IQ( 1,NP )+.5*Q(3,NP 
1) )*RTTO)/((QT(l,NP)+.5*QT(2»NP))-(Q<l,NP)+.5*Q( 2,NP)J*RTT0) )+DINP) 
GO TO 201 

200 CHDLDTI NP)=0. 

ZICTPT(NP)=ZICP(NP) 

ZICTP(NP)=DINP) 

201 CHOLTOI NP) = CHDLD( NP)*CLD£SG/CLT +CHDLDTI NP) 

CLLTOT( NP)=CLL( NP )*CLDESG/CLT+CLB (NP) 

GAMMATt NP)=CHDLDT(NP)/(4.0*QBAR*YMAX) 
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GAMTMA(NP)=GAMMA(NP)*CLDESG/CLT+GAMMAT(NP) 

LOACT { NP )=( ZICTP(NP)— (D ( NP)— C (NP ) ) ) / C 2. 0*C ( NP I) 

LUACTOl NP)=(ZICTPT( NP 1 — t D ( NP )— C ( NP) ) ) / ( 2. 0*C ( NP ) ) 

AMU ( NP ) =— ( ( CHDL D( NP )— P l/2.*(Q(l»NP)+.5*Q(3»NP) ) ) / ( QBAR*4. *C { NP) ) ) 
AMUT ( NP ) =— ( CHDLDT ( NP)-PI/2 . * ( ( QT ( 1 , NP) +.5*QT(3,NP) )-(Q( 1 ,NP )+.5*Q 
1(3,NP) )*RATO)/( QBAR*4.*C(NP ) ) ) 

AMU TOT ( NP ) = AMU ( NP)*CLDESG/CLT + AMU T ( NP ) 

1000 CONTINUE 


AMXB=0. 

ETACP=0 * 

00 715 NP=I SSST »NMAX 
RU=PI*ETA(NP)*VE(NP)/(2.*SSSTA) 
AMXB=AMXB+RU*CHDLDT { NP) 

715 ETACP=£TACP+RU*SLC( NP) 

CMR00T=2.*YREF*AMXB/(QBAR*S«EF) 


CSUCT=0.0 
DO 714 N=1 » NMAX 

RU= YMAX* VE( N)*PI**2/(8.*SSSTA*SREF) 

IF(ETA(N).EQ.O. ) GMICRN=0. 

IF ( ABS ( ETA{ N) ) . GT .0 . . AND. ABS ( ET A ( N ) ) .LE.Y3NP) 0MICRN= I ETAI N ) /ABSl E 
1TA(N) ) )*CHI 

IF( ABS(ETA(N) > .GT.Y3NP.AND. ABS(ETA(N) ).LE.YMIN) OMI CRN= ( ETA ( N ) / ABS 
KETA(N) ) ) ♦ALAME 

IF ( A8S( ETA( N) ) .GT .Y MIN. AND. DELTA. GE .0.) OMICRN= ( ETA (N ) /ABS ( ETA( N) ) 
1 ) *ALAME 

IF( ABS (ETA(N) ).GT. YMIN. AND. DELTA. LT.O.) OMI CRN= (ETA ( N) /ABSCETAJ N) ) 
l)*(PI/2. -DELTA) 

714 CSUCT=CSUCT+RU*BETA*(Q( 1 ,N) /QBAR) **2/ (C (N)*COS( OMICRN) ) 

CDI I=CLT— CSUCT 
489 CCC=0.0 

DO 421 N=1 » NMAX 
CCC=CCC+GAM MA ( N ) **2 
421 CONTINUE 
CCD=0.0 

DO 91 NUP=1 tNMAX 

ALPHAI ( NUP) =GAM MA ( NUP ) / SA ( NUP »NUP ) 

DO 91 N=1 f NMAX 

IF( N.EQ.NUP ) GO TO 91 

ALPHAI (NUP) = ALPHAI ( NUP)-SA( NUP,N) *GAMMA ( N ) / S A ( N UP ,NUP ) 
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NUST=IABSINUP— ISSST I ,■< : . 

IF( MOD( NUST »2)» EQ.O) GO TO 91 
NPST—I ABSIN — ISSST+i| 

IF( MOD( NPST ,2 )• EQ.O ) GO TO 91 
CCD=CCO- 2 -0 *SA( NUP » N i .*4 6AMM A* NUP i *G AMMA C N) l 
91 CONTINUE 

CDI«Pl*CCCC*CCDI/SR£F \ 

DI0L-CDI/CLT**2 

DII0L=CDII/CLT**2 


IF( JTEST.NE.liGO TO 1Q10 
GO TO 1011 

1010 XLREF=XLREF*YMAX1 
CREF=CREF*YMAXi 
REFPT=REFPT*YMAXi 
SREF=SREF*YMAXi**2 

1011 WRITE(6,552) 

WRITE(6,541) ARNN,SREF # REFPT,CREF*XLREF 

WRITE OUTPUT DATA 

WRI TE (6 *510 ) 

WRI TE(6»502 ) 

WRITE (6,530) CASE 

WRITE (6, 559) ALPHZO»CMZ£RO,CMROOT 

WRITE (6,546) CLT,CH,CRL ,CLTWST,CMTAC,CLOESG 

WRI TE (6,531 ) CLTTfCMTfD IOL ,OIIOL 

WRI TE(6 ,554 ) 

WRITE (6, 52 8) ZICPT,TOAC,ETACP 
WRI TE(6 ,510 ) 

WRITE(6*520) 

WRI TE(6 ,521 ) 

WRI TE (6 ,538 ) 

WRITE (6, 547) 

DO 850 N=ISSST,NMAX . 

850 WRI TE (6 ,511 ) N,ETA(N) , GAMMA! N) ,GAMMATt N) , GAMTMAI N) , CLL (N ) , CLB( N) ,CL 
1LT0T (N) , AMU (N) , AMUT ( N) , AMUTOT (N) ,C( N) ,D( N) 

WRI TE (6 ,510 ) 

WRITE (6,542) 

C 

DO 188 N=ISSST , NMAX 
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DO 188 LSv.lf.lD 
ALS=LS 

TAW(LS) =ALS#PI/ 10.0 
ZI=(1.— COSITAW(LS)) )/2. 

TAT=0.0 

TATWST=0. 

DO 185 J=2» JMAX 
A J= J 

TAT = TAT+Q(J»N)*SIN( ( AJ-i .0 ) *T AW (L S ) ) 

TATWST=TATWST+( QT ( J , N)-Q( J , N) *R ATO) *S IN H AJ- 1 .0 ) *TA W( LS ) ) 

185 CONTINUE 

PAU=C0S(TAW(LS)/2.G)/SIN(TAW{LS)/2.0) 

DEL TP=1./C(N)*(Q(1,N)*PAU+T AD 

DELPTW=1./C<N)*((QT(1,N)-Qi 1 ,N) *RATO) *P AU+TATWST) 

DLPTWC=DELT P*CLDESG/CLT +DELPTW 

WRI TE ( 6 , 543 ) N, LS ,E TA ( N ) , Z l , OELTP ,DELPT W, OLPTWC 
188 CONTINUE 
C 

WRI TE(6 ,510 ) 

DO 1060 K=i y JMAX 
1060 J J J ( K ) = K— 1 

WRITE(6,516 ) 

WRI TE(6,519) 

WRI TE( 6,539 )( JJJ< I) ,1=1, JMAX) 

DO 1030 N=I SSST , NMAX 

1030 WRI TE l 6 , 534 ) N, ( Q II , N ) , I =1 , JMAX ) 

WRI TE ( 6 , 523 ) 

WRITE (6,519) 

WRITE(6,539)(JJJ( I) ,1=1, JMAX) 

DO 1031 N=I SSST, NMAX 

1031 WRI TE ( 6 , 534 ) N, ( QT( I , N) , 1 = 1 , JMAX ) 

383 WRI TE(6,510) 

WRI TE ( 6 ,517 ) 

WRITE(6,518) 

WRI TE(6,512) 

WRI TE ( 6 , 513 ) 

DO 378 N=ISSST , NMAX 

378 WRI TE(6,533) N, ETA( N) ,CHDLD (N) , CHDLOT IN ) ,CHDLTO(N) , ZICP ( N ) , ZI CTP ( N 
1 ) . 7 1 CTPT ( N) .LOAC(N ) ,LQACJ IN ) ,LOACTO( N ) , SLC( N) , BSPLDTIN ) ,8LDT( N) 


GO TO 1 
10 STOP 
END 



FUNCTION FOFT1 ( THET A) 

COMMON XSUB » YSUB » ET ASUB , CSUB , DSUB 
XI=—CSU8*C0S ( THETA) +OSUB 

3K= I. +( XSUB-XI) /SQRTI IXSUB-XI ) **2 + ( YSUB-ETASUB ) **2) 
FOF T1-BK*( 1 «0+C OS (THETA ) ) 

RETURN 

END 


FUNCTION FOFT2I THETA) 

COMMON XSUB , YSUB »ET ASUB ,C SUB, DSUB 
XI=— CSUB*COS( THETA) +OSUB 

8K=L.-M XSUB-XI) /SQRTI (XSUB-XI ) **2 +1 YSUB-ETASUB) **2) 

FOFT2— BK*SI N( THETA) **2 

RETURN 

END 


FUNCTION F0FT3I THETA) 

COMMON XSUB , YSUB, ETASUB, CSUB, DSUB 
XI=-CSUB*COS( THETA) +DSUB 

BK=I . + { XSUB— X I ) /SQR T ( (XSUB-XI ) **2+I YSUB-ETASUB) **2) 
FOFT3=BK*SI N( THETA) *S IN ( 2 .0 *THET A) 

RETURN 

END 


FUNCTION F0FT4I THETA) 

COMMON XSUB, YSUB, ETASUB, CSUB, DSUB 
XI=-CSUB*COS( THETA) +DSUB 

BK=1 . + ( XSUB-Xl )/SQRT( ( X SUB- X I ) **2+< YSUB-ETASUB ) **2 ) 
F0FT4=BK*SI Nl THETA) *S IN I 3 .0 *THET A) 

RETURN 

END 


FUNCTION F0FI5I THETA) 

COMMON XSUB, YSUB, ETASUB , CSUB, DSUB 
Xl=-CSUB*CQSl THETA) +DSUB 

BK-l.+I XSUB— XI ) /SORT l I XSUB-X I) **2 + l YSUB-ETASUB ) **2 ) 
F0FT5=BK*SIN(THETA)*SIN{4.0*TH£TA) 

RETURN 

END 



FUNCTION F0FT6I THETA) 

COMMON XSUB , YSUB ,ET ASUB ,C SUB, DSUB 
XI =— CSUB*CO S ( THET A) +DSUB 

BK=l.+( XSUB-XI ) /SORT ( (XSUB-XI )**2+( YSUB-ETASUB ) **2) 
F0FT6=BK*SI N( THETA ) *S IN ( 5 . 0 *TH£ TA ) 

RETURN 

END 


FUNCTION FOFT71 THETA) 

COMMON XSUB , YSUB , ET ASUB ,CSUB , DSUB 
XI=— CSU8*C0S( THETA) +DSUB 

BK=l.+( XSU8-XI ) /SQRT( (XSUB-XI )**2+( YSUB-ETASUB) **2) 
F0FT7=BK*SI N( THETA) *SIN(6»0*THETA) 

RETURN 

END 


FUNCTION FOFT8I THETA) 

COMMON X SUB, YSUB,ETASUB,CSUB, DSUB 
XI CSUB*COS( THETA) +DSUB 

BK= 1 . + ( XSUB-XI ) /SORT ( (XSUB-XI )**2+( YSUB-ETASUB ) **2) 
FOF T 8=BK*SI N( THETA) *S IN { 7 .0 *THETA ) 

RETURN 

END 


FUNCTION F0FT9( THETA) 

COMMON XSUB*YSUB,ETASUB,C SUB, DSUB 
XI=— CSUB*COS( THETA) +DSUB 

BK= I . + ( XSUB-XI ) /SQR T ( (XSUB-XI )** 2+1 YSUB— ETA SUB ) **2) 
F0FT9=BK*SIN(THETA)*SIN(8.0*THETA) 

RETURN 

END 


FUNCTION FOFT 10 (THETA ) 

COMMON XSUB, YSUB, ET ASUB »C SUB, DSUB 
XI=-CSUS*COS{ THETA) +DSUB 

BK=1 .+( XSUB-XI ) /SORT ( (XSUB^Xl )**2+( YSUB-ETASUB) **2) 
F0FT10=BK*S IN(THETA)*SIN(9.0*THETA) 

RETURN 

END 



o o o o no 


SUBROUTINE GAUSS I A , B,N , SUM ,FQFX ) 

REFERENCE SCARBOROUGH NUM. MATH. ANAL. PAGE 147 
HOWEVER THIS SUBROUTINE INTEGRATES FROM ZERO TO ONE 

DIMENSION U ( 5 ) » R (5 ) 

U( I >. = .425562830509184 
U< 2 >=.283302302935376 
U( 3 ) = . 16029 5215 850488 
UI4)=. 06746 83 16 655508 
U( 5 )=.0 13046735741414 
R{ L )=. 14776 2 L 12 35 7376 
R C 2 ) — • 1 3463 3359654998 
R< 3) = . 109543181257991 
R (4 >=.074725674575290 
R(5 >=.033335672 154344 
SUM=0.0 

IF(A.EQ.B) RETURN 
F I NE=N 

DEL TA=F INE/ ( B-A ) 

DU 3 K=1,N 
XI=K— 1 

F INE = A+XI /DELTA 
DO 2 11= 1,5 
UU=U(I I J/DELTA+FINE 

2 SUM=R(II) *FOFX(UU)+SUM 
DO 3 L = 1 ,5 

UU=( l.O-U(L) ) /DELTA+F INE 

3 SUM=R(L)*FOFX(UU>+SUM 
SUM= SUM/DELTA 
RETURN 

END 

MATRIX INVERSION WITH ACCOMPANYING SOLUTION OF LINEAR EQUATIONS 
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o n o o o o o o o 


c 


SUBROUTINE MATI NV ( A , N, B , M, DETERM, I P I VOT , INDEX, NMAX, ISCALE) 


DIMENSION I PI VOT (N) , A (NMAX, N) ,B INMAX, M) , INDEX! NMAX, 2) 
EQUIVALENCE ( IROW,JROW) , { ICOLUM, JCOLUM), (AMAX, T, SWAP) 

INITIALIZAT ION 

5 I SCALE=0 

6 Rl=10. 0**13 

7 R2=l. 0/R1 
10-DETERM=1.0 
15 DO 20 J = 1,N 
20 I PI VOT ( J ) -0 
30 DO 550 1=1, N 

SEARCH FOR PIVOT ELEMENT 

AO AMAX=0. 0 
A 5 DO 105 J=1 , N 

50 IF ( IPI VD1 ( J)-l ) 60, 105, 60 
60 DO 100 K=l,N 

70 IF ( IPI VOT(K)-i) 80, 100, 7A0 

80 IF ( ABS(AMAX)-ABS(A(J,K) ))85,100,100 

85 IROW=J 

90 ICOLUM=K 

95 AMAX=A(J,K) 

100 CONTINUE 

105 CONTINUE 

IF (AMAX) 110,106,110 

106 DETERM=0»0 
I SCAL£=0 
GO TO 7A0 

110 IPI VOT ( ICQLUM)= IPIVOT ( ICOLUM) +1 

INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL 

130 IF ( IROW-ICOLUM) 1A0, 260, 1A0 
1A0 DETERM=-DET ERM 
150 DO 200 L=1 , N 
160 SWAP=A( IROW ,L ) 

170 A( I ROW, L )=A ( ICOLUM, L) 

200 A( ICOLUM, L)=SWAP 
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non non 


205 IF ( M ) 260, 260, 210 
210 DO 250 1=1, M 
220 SWAP=B< IR0W,L) 

2 30 6! I ROW, L ) =B ( I COLUM, L ) 

250 8! ICOLUM, L)=SWAP 
260 INDEX! I , 1 )= IROW 
270 INDEX! I »2) = ICOLUM 
310 PIV0T=Al ICOLUM, ICOLUM) 

SCALE THE DETERMINANT 

1000 PIVOT I=P 1 VO T 

1005 IF! ABS! DETE RM )-Rl ) 1 030, 1010, 1010 
1010 0ETERM=DETERM/R1 
ISCALE=ISCALE+1 

IF! ABS! DETERM )-Rl) 1060, 1020,1020 
1020 DETERM=DET£RM/Rl 
I SCALE= I SCALE + 1 
GO TO 1060 

1030 IF! ABS! DETERM )-R2) 1040, 1040,1060 
1040 DE TER M= DETERM #R1 
ISCALE=ISCALE— 1 

IF! ABS! DETERM) — R2)l 050, 1050,1060 
1050 DETERM=DETERM*R1 
I SC ALE= I SCALE-1 

1060 IF! ABS! PIVOT I )-R 1)1 090, 1070,1070 
1070 PIVOT I=PI VO TI/R1 
ISCALE= fSCALE+1 

I F ! ABS! P I VOT I >-R 1)320, 1080, 1080 
1080 PIV0T1=PIV0TI/R1 
ISCALEMSCALE+1 
GO TO 320 

1090 IF! ABS! PIVUTI )-R2 )2000, 2000 , 320 
2000 PIV0TI=PIV0TI*R1 
ISCALE= ISCALE— 1 

IF! ABS! PI VOT I )-R2 )20 10, 2010,320 
2010 P I V 0 1 1 = P I VO T I *R 1 
I SC ALE= I SCALE— 1 
320 DETERM=D£TERM*P IVOTI 

DIVIDE PIVOT ROW BY PIVOT ELEMENT 

330 A! ICOLUM, ICOLUM)=1»0 



non ocvo 


340 DO 350 L=l,N 

350 A( IC0LUM,L)=A( 1C0LUM, D/PIVOT 
355 IF ( M ) 380, 380, 360 
360 DQ 370 L=l, M 

370 BUC0LUM,L} =B( IC0LUM,LJ /PIVOT 

REDUCE NON- PIVOT ROWS 

380 DO 550 L1=1,N 

390 IF( Ll-ICOLUM) 400, 550, 400 

400 T=AlLI, ICOLUM) 

420 Al L 1 » ICOLUM )=0» 0 
430 DO 450 L=I , N 

450 A(LI,L)=A(L I ,L>— At I COLUM, L) *T 
455 IF ( M ) 550, 550, 460 
460 DO 500 L-l » M 

500 B(L1,L)=B(LI,L)— B(I COLUM, L ) *T 
550 CONTINUE 

INTERCHANGE COLUMNS 

600 DO 710 1=1, N 
610 L=N+1-I 

620 IF ( INDEX! L , 1 )-INDEX(L, 2) ) 630, 710, 630 
630 JROW=INDEX( L , 1 ) 

640 JCO L UM= I NDE X ( L , 2 ) 

650 DO 705 K=1,N 
660 SWAP=A( K,, JROW) 

670 A( K , JRQ W ) = A (K » JCGLUM) 

700 AIK , JCOLUM) =SWAP 
705 CONTINUE 
710 CONTINUE 
740 RETURN 
END 
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Sample Output Listing 
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HAS A TRAILING EDGE CHORD EXTENSIONI 
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CASE NUMBER- 5 

ANGLE FOR ZERO LIFT* ALPHA ZERO* -0.00000 PITCHING MOMENT CQEFF. AT ZERO LIFT*CM ZERO* 0.00000 

ROOT BENDING MOMENT COEFFICIENT AT ZERO LIFT»CBMROOT* 0.00000 
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MEAN CAMBER SURFACE PROGRAM A0457 
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Program Listing 



PROGRAM MNCAMBR I INPUT , OUTPUT ,TAPE5= INPUT, TAPE6=0UTPUT ) 

TO SOLVE FOR THE MEAN CAMBER SURFACE REQUIRED TO SUPPORT A GIVEN 
PRESSURE LOADING ON A VARIABLE SWEEP PLANFORM, ( Z/C ) VS. (X/C) 


DIMENSION Y{41) ,ETA(41), CI4l),D( 41) ,QP( 110* li ,PHlUO) , VEI41) , B( 41) 
UAL ( 10,41) ,DF( 101 ,SUML( 110, 110 ) ,YUT { 10,41) , W( 10 ) ,RESID( 10,41) , SUMS 
2T(41) , AST (10,10 )»CUT( 10,10) , CONST ( 1 10 , 1 ) 

REAL MACH 

COMMON XSUB » YSU8 »ETASUB ,CSUB»DSUB ,B SI 10, 1 ) , X( 10) ,XC,NS 
EXTERNAL F0FT1 , F0FT2 , FOFT 3, FOFT 4 , F0FT5, F0FT6 ,F0FT7» F0FT8,F0FT9,F0F 
ITIO , WU 

500 FORMAT (6F12..5) 

501 FORMAT l 1HI , 58 X, 13HGE0METRY DATA///) 

502 FORMAT ( 1H1»57X» 16HMEAN CAMBER DATA///) 

503 FORMAT ( 50X, 5HX2 PP=, F9 .5 , 5X , 5HY2 PP=, F9 .5 ) 

504 FORMAT ( 50X, 5HX3 NP= , F9. 5 , 5 X, 5HY3NP=, F9.5 ) 

505 F0RMAT{50X,5HX4PP=»F9.5»5X,5HY4PP=»F9.5 ) 

506 FORMAT! 50X, 5HX5PP=,F9. 5 ,5X, 5HY5PP=, F9.5 ) 

507 FORMAT ( 50X, 5 HX6 NP= ,F9.5,5X,5HY6NP=,F9.5) 

508 FORMAT! 50X, 5HX7 PP= , F9. 5 ,5X, 5HY7PP=,F9.5 ) 

509 FORMAT !40X1 2HCA SE NUMBER=*F6.0,5X*14H SYMMETRY CODE* ,F5.0,5X, 12HMAC 
1H NUMBER= ,F 9. 5 ) 

510 FORMAT! 1H0) 

511 FORMAT ( 8F9. 5) 

512 FORMAT! 1H1) 

513 F0RMAT(6X3HX/C9X3HX/C9X3HX/C9X3HX/C9X3H X/C9X3HX/C9X3HX/C9X3HX/C9X3 
1HX/C9X3HX/C ) 

514 FORMAT! 15X,97HIF SYMMETRY CODE IS EQUAL TO l, THE SPAN LOAOING IS S 
1YMMETRI CAL/ OTHER THAN 1 , IT IS ANT I S YMME TR IC AL ) 

515 FORMAT! ////52X, 17HMEAN CAMBER SHAPE//) 

516 FURMAT(4XF7.5»5XF7.5»5XF7.5»5XF7.5»5XF7.5»5XF7.5f5XF7.5»5XF7.5»5XF 
17.5 ,5XF7.5) 

517 FORMAT! 55 X, 3HX/C , 6X , 3HZ/C ) 

518 FORMAT ( 20X, 15HSTATI ON NUMBER= , 1 4 , 10X, 18 HSPANW I SE LOCAT 1 0N=, F9 . 5 , 10 
1 X , 6HCH0RD= , F9. 5 ) 

519 FORMAT! 10F12.5) 

520 FORMAT! 7X,2HA1 , 1 OX , 2HA2 , 10X , 2H A3 , 10X ,2HA4, 10X , 2HA5 , 10X , 2HA6, 10X2H 
1 A7 , 10X, 2HA8 , 10X , 2HA9 , 10X , 3HA 1 0) 

521 FORMAT! 4X, 8HC0N.PT. 1,4X, 8HC0N.PT. 2, 4X , 8HC0N.PT. 3,4X ,8HC0N.PT.4, 4X, 
18HC0N.PT.5, 4X, 8HC0N.PT. 6, 4X, 8HC0N.PT. 7, 4X, 8HCQN.PT. 8, 4X,8HC 0N.PT.9 
2, 4X, 9HC0N.PT. 10) 
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522 FORMAT l 30X1 3HASPECT R AT 1 0=, F9. 5 * 5X, 14HPLANF GRM ARE A=, F9 . 5 , 5 X, 14HAV 
1ERAGE CHORD = f F9 .5//10X21HMEAN GEOMETRIC CHORD= , F9 . 5 , 5X, 59HX LOCATI 
2 ON OF THE LEADING EDGE OF THE MEAN GEOMETRIC CH0RD=,F9.5/ 30X39HY 
3L0C AT ION OF THE MEAN GEOMETRIC CHORD= , F9. 5// ) 

524 FORMAT! 15X, 1IHR00T CH0RD=,F9.5»5X,l OHTI P CHORD= , F9. 5, 5X , I5HF0RE W IN 
1G CH0RD=,F9.5,5X,20H0VERALL TAPER RATIQ-»F9.5//45X11HY LE BREAK=, F 
29.5 , 5X, 1 IHY TE BREAK= »F9 . 5/ /25X17HX PIVOT LOCAT ION= ,F9. 5 , 5 X17HY PI 
3 VOT L0CATI0N=,F9.5,5X,17HZ PIVOT LOCATION^ F9 . 5// 55 X, 19HTE CHORD E 
4XTE NS ION= ,F 9 . 5 ) 

525 FORMAT ( 4F6. 0) 

526 FORMAT! ////40X5 9HL0CATI0N OF PERIMETER POINTS FOR- THE PLANFORM USE 
ID A§ INPUT) 

527 FORMAT! //5X35HNUMBER OF CHORDWISE PRESSURE MODES= »F5.0, 5X72HNUMBER 

1 OF STATIONS SPANWISE ON A PANEL WHERE PRESSURE MODES ARE DEFINED^ 
2,F5.0) 

528 FORMAT! 40X, 49HSL0PE S» (W/U) »AT CONTROL POINTS, FROM FRONT TO REAR//) 

529 FORMAT t 50X, 9HX1PP= 0.10X,9HYIPP= 0.) 

530 FORMAT ! 57X, 12HCASE NUMBERS » F6. 0 ) 

532 FORMAT! 52X»7HXi= 0.,10X,7HYl= 0.) 

535 FORMAT ('52 X, 3HX2=,F9.5, 5X, 3HY2=,F9.5 ) 

536 FORMAT! 52X, 3HX3= , F9 . 5 , 5 X, 3H Y3= , F9 .5 ) 

537 FORMAT! 52 X, 3HX4= , F9 . 5 , 5 X , 3HY4= , F9 .5 ) 

540 FORMAT ( 50Xt 5HX6AP=, F9.5 , 5X, 5HY6 AP=, F9.5 ) 

545 F0RMAT(52X,3HX5=,F9.5,5X,3HY5=,F9.5 ) 

548 FORMAT! 5X,23HLE INBOARD SWEEP ANGLE= , F9 . 5 , 5X , 32 HL E INITIAL OUTBOAR 
ID SWEEP ANGL£=, F9.5 ,5X» 30HLE FINAL OUTBOARD SWEEP ANGLE=»F9.5/5X f 
223HTE INBOARD SWEEP ANGLE=» F9.5*5X, 32HTE INITIAL OUTBOARD SWEEP AN 
3GL£=»F9.5»5X| 30HTE FINAL OUTBOARD SWEEP ANGLE= » F9 .5//5X ,40 
4HCHANGE IN OUTER PANEL SWEEP ANGLE , DELT A= , F9. 5 , 4X , 26HP I VOT CANT AN 
5GLE IN PITCH=,F3.0,4X,25HPI VOT CANT ANGLE IN R0LL=,F3.0// ) 

549 FORMAT ( 5 IX, 4HX6 A= ,F9. 5 , 4X ,4HY6A= , F9 . 5 ) 

550 FORMAT! 52 X, 3HX6= , F9 . 5 , 5 X , 3HY6= , F9 . 5 ) 

551 FORMAT! 2X,122HTOTAL WING PLANFORM (MEAN GEOMETRIC CHORD AND ITS LOC 
I AT I ON ARE INVALID IF THE PLANFORM HAS A TRAILING EDGE CHORD EXTENS 

2 ION)///) 

553 FORMAT! ///3 3X ,7 2HL0CAT I ON OF PERIMETER POINTS FOR PLANFORM TO BE U 
ISED IN THE COMPUTAT I ONS /40X , 57HWHEN NONOIMENSIONALI ZED BY THE SEMI 
2 SPAN RATIO GIVEN ABOVE) 

555 FORMAT ( 40 X, 57HWHERE THE ORIGIN IS AT THE HALF ROOT CHORD, POSITI VE 
IX AFT//) 

556 FORMAT! 52X, 3HX7=, F9 .5 , 5X, 3H Y7= , F9.5 ) 

558 FORMAT! 25X, 33H! SEMI SPAN AT FINAL OUTBOARD SWEEP, 37H/SEMISPAN AT IN 
1ITIAL OUTBOARD SWEEP )=, F9. 5 ) 
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n n n n n n n n <r. <n c) n n n n n n c> n n n n n ri 


559 FORMAT ( 31 X6 9HP0L YNOMI AL COEFFICIENTS OF (X/C),IN OROER OF INCREASI 
1NG POWERS OF X/C//) 

560 FORMAT! 51X,2F9. 5) 

INPUT DATA 


IF THE WING DOES NOT HAVE A LEADING EDGE BREAK SET BIRAT AND CHI 
EQUAL TO ZERO 


IF THE WING DOES NOT HAVE TRAILING EDGE BREAK SET B2RAT AND PSI 
EQUAL TO ZERO 


IF THE WING DOES NOT HAVE A VARIABLE SWEEP OUTER PANEL SET DELT, 
XP AND VP EQUAL TO ZERO 


IF THE SPAN LOADING IS TO BE SYMMETRICAL SET THE SYM CODE EQUAL TO 
ONE, IF ANTI SYMMETRICAL SET THE SYM CODE EQUAL TO TWO 


I RE A D ( 5 » 5 0 0 ) AR , CHI , AL AMD , PS I , BIRAT, B2RAT 
IF! EOF, 5) 3,4 

4 READI5, 500) TAP ER , DELT, XP ,YP,CHDEXT , MACH 
READ! 5,525) CAS E , SYM,CST A, SSST A 


VARIBLE SWEEP GEOMETRY PROGRAM 


PI=3. 14159265 
RAD=1 80 . /P I 
QBAR=1. 00 
B0= 1.0 
YMI N= 1.0 
YMAX=1 . 0 
YMAXl = l .0 
SIGM=0. 

RHH=0. 
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c 


c 

c 


ZP=0. 


CHDEXX=CHDE XT 
BETA=SQRT ( I MACH**2) 

I SYM=SYM 

JMA X=CST A 

MS=CSTA 

NS=CSTA 

I SSST=SSSTA 

NMAX=2* I SSS T-l 

NNI I=ISSST+ 1 

JKMAX= J MAX* ISSST 

MSMAX=10 

NSMAX=10 

LS= 1 

ITTU= 1 

CHI =CHI /RAD 
AL AM D= ALAMO /RAD 
PSI =PSI /RAO 
DELTA=DELT/RAD 
SIGMA=S IGM/RAO 
RHO=RHH/RAO 
TANC=TAN(CHl ) 

TANL=TAN( ALAMO) 

TANP=TAN (PS I ) 

TANDE=T AN( P 1/2. +OEL TA ) 

B1=B1RAT*B0 

B2=B2RAT*BO 

CR=BO*( 4./AR-B2RAT*TANP-TANC* (BIRAT *( BIRAT-B2RAT-1. ) )-T ANL* 

1 ( BI RAT* ( B2R AT-BIRAT+i .)-B2RAT ) ) *( 1. /( B2RAT* ( 1. -TAPER )+( l.+TAPER) ) ) 
OMEGA=ATAN ( ( 1 . / ( 1. -B2R AT ) ) *( I TAP ER-1 . ) *( CR/BO ) +B 1RAT* (TANC-TANL ) 
L-B2RAT*TANP+TANL ) ) 

TANO=TAN( OMEGA) 

EOMEG=OMEGA+DELTA 
ALAME=ALAMD+OELTA 
TANE=TAN ( EOMEG) 

TANA=TAN( ALAME) 

CR=CR+CHDEXT 
X2=— CR/2. 

Y2=0.00 
Z2=0. 00 
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X3=-CR/2.+Bl*TANC 

Y3=B1 

Z3=0.00 

X4=-CR/2.+Bl*(TANC-TANL)+B0*TANL 

Y4= BO 

7.4=0.00 

X5=CR/2 . J-62*I TANP— T ANO)+BO*TANO— CHDEXT 
Y5= BO 
Z5=0 . 00 

X6A=CR/2. +82+TANP 

Y6A=82 

Z6A=0.00 

X6=X6A-CHDEXT 

Y6= Y6A 

Z6= 0.00 

X7=CR/2. 

Y7=0.00 

Z 7=0. 00 

X2P P=X2 

Y2PP=Y2 

Z2PP=Z2 

X6AP=X6 A 

Y6AP=Y6A 

Z6AP=Z6A 

X7P P=X7 

Y7PP=Y7 

Z7PP=Z7 

IF! CELT A.EQ.O. ) GO TO 360 

X3PP = XP+IX3-XP)*COS( SIGMA )*COS I DELTA )+( Y3-YP )*COS ( RHO)*SIN( DELTA) 
1 +( Z3-ZP ) * ICQS I S IGMA )*SI NI RHO ) *S INIDELTA )— SIN I SIGMA ) *COS( RHO ) ) 
Y3PP=YP — (X3 — XP) *COS( SIGMA )*S INI DELTA )+( Y3-YP ) *COS I RHO)*COS( DELTA) 
1 + 1 Z3-ZP ) *( COS I S IGMA )* SI N (RHO) *GOS( DELTA ) +SINI S IGMA )*SIN( DELTA) ) 
Z3PP= I X 3-XP ) * SI N ( S I GMA ) *COS I RHO ) — I Y 3— YP ) *S I N( RHO ) 

1+IZ3-ZP )*COSISIGMA)*COS(RHO)+ZP 

X4P P = XP+ ( X4— XP) *COS(SIGMA)*COS(DELTA)+( Y4-YP )*COS I RHO )*S IN I DELTA) 
1+ 1 Z4-ZP ) *(COS ( S IGMA )*S INI RHO ) *S INI DELTA ) — S IN IS IGMA )*C0SIRH0)) 
Y4PP=YP — (X4— XP) *COS(SIGMA)*SIN(DELT A) +( Y4— YP ) *COS I RHO)*COS( DELTA) 
I + IZ 4-ZP ) * (COS ( SIGMA ) * S I N I RHO ) *COS( DELTA) +SINI S IGMA ) *SI N( DELTA ) ) 
Z4P P= I X4- XP ) * S I N I S I GMA ) *COS I RHO ) - ( Y4- YP ) * S I N I RHO ) 

1+ 1 Z 4-ZP ) *C0 S I S I GMA ) *COS I RHO ) +ZP 

X5PP=XP + (X5 — XP)*COS( SIGMA)*COS(DELTA)+( Y5-YP)*C0S(RH0)*S INI DELTA) 
1 + (Z5-ZP)*ICC)S(SI GMA ) *S I N I RHO )*SIN(DELTA) — SINISIGMA) *COS( RHO) ) 
Y5PP=YP— ( X5— XP) *COS( SIGMA)*SIN(DELTA)+( Y5-YP )*COSI RHO)*COSI DELTA) 
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1-MZ5-ZP)*(CQS(S IGMA ) * S I N ( RHO ) *COS ( OELT A) +S IN ( S IGMA ) *S IN ( DELTA ) ) 
Z5PP= { X5-XP ) *SIN( SIGMA ) *CQS ( RHO ) - ( Y5- YP ) *51 N I RHO ) 

1+ (Z5-ZP)*COS(S IGMA) *COS (RHO )+ZP 

X6PP=XP+ ( X6— XP) *C0S ( SIGMA) *C0S( OELT A) +( Y6-YP)*COS4RHO)*SINCOELTAJ 
1+IZ6-ZP)*(CDS(S IGMA ) *5IN( RHO) *S INC DEt TAJ- S I N { S I GM A ) *G OS (RHO J ) 
Y6PP=YP— (X6— XP) *COS( SIGMA )*SIN(DELT A)+(Y6-YPJ*C0S I RHO )*COS( DELTA) 
I+(Z6-ZP)*(C0S( SIGMA)* SIN(RHO)*C0SI DELT A ) +SIN( S IGMA )*SIN( DELTA) ) 
Z6P P= ( X6-XP ) * S I N ( SI GMA ) *COS (RHO )- { Y6- YP ) *S I N ( RHO ) 

I + ( Z6-ZP ) *COS( SIGMA) *COS (RHO J+ZP 
IF ( ( Y3-Y2PP ) . EQ .0. ) GO TO 208 
A2PP=(X3 -X2PP ) / ( Y3 -Y2PP) 

GO TO 209 

208 A2PP=0. 

209 IP C ( Y6-Y7PP ).EQ.O.) GO TO 210 

A6PP = (X6 - X7PP ) / ( Y6 -Y7PP ) 

GO TO 211 

210 A6PP=0. 

211 IF( ( Y4PP-Y3 PP) . EQ.O . ) GO TO 212 
A3PP= ('X4PP-X3PP ) / (Y4PP-Y3PP ) 

GO TO 213 

212 A3 PP=0. 

213 IF( (Y5PP-Y6PP).EQ.0.) GO TO 214 
A5PP= ( X5PP-X6PP ) / (Y5PP-Y6PP ) 

GO TO 215 

214 A5PP=0. 

215 IF ( ( { X3— X2P P ) ** 2 + ( Y3— Y2PP ) **2 ) » EQ .0 . ) GO TO 216 
G2PP=(Z3 - Z2PP ) / ( ( X3 — X2PP)**2+(Y3 - Y2PP ) **2 ) **. 5 
GO TO 217 

216 G2PP=0. 

217 IF ( ( ( X7PP-X6) **2 + (Y7PP-Y6 )**2 ) . EQ.O » ) GO TO 218 
G6PP=(Z6 -Z7PP ) / ( ( X7PP-X6 ) **2+ ( Y 7PP- Y6 )**2)**.5 
GO TO 219 

218 G6P P=0» 

219 IF ( ( A3PP— A2 PP ) . EQ.O . ) GO TO 220 

X3NP= ( A3PP*X2PP-A2PP*X4PP+A2PP*A3PP*(Y4PP-Y2PP ) ) / ( A3PP-A2PP ) 
Y3NP= (X2PP-X4PP+A3PP*Y4PP-A2PP*Y2PP)/( A3PP-A2PP) 

GO TO 221 

220 X3NP=X2PP 
Y3NP=Y2PP 

221 Z3NP=< ( X3NP— X2PP )** 2+( Y3NP— Y2PP ) **2 ) **. 5*G2PP+Z2PP 
IF( ( A6P P— A5 PP ) • EQ.O • ) GO TO 222 

X6NP= ( A6PP*X5P P-A5PP*X7PP+A5PP*A6PP*( Y7PP-Y5PP ) ) / ( A6PP-A5PP ) 
Y6NP=( X5PP-X7PP + A6P P*Y7PP— A5PP*Y5PP ) / ( A6PP-A5PP ) 
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GO TO 223 

2 22 X6NP=X7PP 
Y6NP=Y7PP 

223 Z6N P= l ( X6NP-X7PP )**2MY6NP”Y7PP )*^2)*** 5*G6PP+Z7PP 
IFCCHOEXT.eg.O. ) GO TO 342 

X6NP-X6PP+C Y6-Y6PP} *A5PP 
Y6N P=Y6 

Z6NP= < X6AP-X6NP > * ( Z6A-Z6 ) / C X6A-X6 J + Z6AP 

342 IF ( Y4PP *GE* Y5PP ) GO TO 411 
IF { Y4PP.LT. Y5PPI YM AX-Y5PP 
YMAX 1=YM AX 
GO TO 4 LQ 

411 YMA X=Y4 PP 
YMAX1—Y MAX 

410 X2PP=X2 PP/YMAX 
Y2P p- Y2 PP/YMAX 
Z2P P = Z2 P P/Y MAX 
X 3N P = X 3 NP / Y M A X 
Y3NP=Y3NP/YMAX 
Z3NP=Z3NP/YMAX 
X4PP=X4PP/Y MAX 
Z4PP-Z4PP/Y MAX 
X5P P=X5 PP/YMAX 
Z5P P=Z5 PP/Y MAX 
X6A P = X6 AP /Y MAX 
Y6AP=Y6 AP/Y MAX 
Z6AP=Z6AP/YMAX 
X6N P= X6 NP /Y MAX 
Y6NP—Y6 NP/Y MAX 
Z 6 N P= Z6 NP / Y MA X 
X7PP=X7PP/Y MAX 
Y7PP=Y7 PP/YMAX 
Z7PP=Z7 PP/YMAX 
XP= X P/Y M AX 
YP=YP/YMAX 
ZP= Z P/Y MAX 
CR= CR / Y MAX 
Y4PP=Y4 PP/YMAX 
Y5P P= Y5 PP/YMAX 
CHDEXT=CHOe XT/YMAX 

IF 1 ABS1 AL AME-CHI ) *G£*0* 000174 ) GO TO 765 
X3N P-XZPP 
Y3N P = Y2 PP 
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Z3NP-Z2PP 

765 IF ( A8S( EOMEG— PS I ) •GE»0« OOOl 74.0R* CHQE XT »NE» O.O) GO TO 767 
X6NP-X7 PP 

Y6N P=Y7 PP 
Z6NP=Z7PP 
X6AP=X7PP 
Y6A P=Y7 PP 
Z6AP=Z7PP 

767 IF ( CHOEXT*NE .0*0) GO TO 766 
X6AP=X7PP 
Y6AP-Y7PP 
Z6A P=Z7 PP 

766 IF C Y6AP * EQ* Y6NP ) GO TO 1050 
CHDEXX=0.0 

GO TO 1051 

1050 CHDEXX-X6AP— X6NP 

1051 SADD=CHDEXX*Y6NP 
YMAX=1.00 

YMI N- AM INI ( Y4PP , Y5PP ) 

81R AP-Y3NP 
82R AP-Y6NP 

IF! Y6NP-Y7PP) 362,361,361 
362 X6NP1=X6NP 
Y6NP1-Y6NP 
Z6NP1—Z6NP 

X6NP-X6NP1- Y6NP1*TANE 
Y6NP=0. 0000000 

T AM U = I Z 6 N P L - Z 5 P P ) / U X 6 N P 1 - X 5 P P ) .2+ ( YoN P 1- Y 5 P P ) ** 2 ) ** , 5 

Z6NP= Z6NPI-I ( X6NP-X6NP 1 ) **2 + ( Y6NP-Y6NP 1 ) **2 ) ** .5*T ANU 

CR=X6NP-X2PP 

ORI GNN- ( X6NP—X7PP ) / 2« 

X2PP=X2PP-ORIGMN 
X3NP-X3NP-ORIGNN 
X4PP=X4PP-0RIGNN 
X5PP-X5.PP-0RI.GNN 
X6NP=X6NP-0RIGNN 
X6AP=X6AP~0RIGNN 
X7PP=0« 00 
B2R AP=0 * 00 
TANP=0. 000000 
GO TO 361 
360 X3NP=X3 
Y3NP-Y3 
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Z3NP=Z3 
X4PP=X4 
Y4P P-Y4 
Z4P P=Z4 
X5PP=X5 
Y5PP=Y5 
Z5PP=Z5 
X6NP=X6 
Y6NP-Y6 
Z6NP=Z6 

SA0D=CHDEXX*Y6NP 
BIRAP=Y 3NP 
B2RAP=Y6NP 
361 TA.NOP=TANE-TANP 
TANUL=T AN E- TANA 
TANOC=TAN£-TANG 
TANLC=TANA- TANC 
TANLP=T ANA-TANP 
TAN PC = T A NP- T A NC 
IF C Y4PP.GE.Y5PP ) XC8LEM= 

1(CR+Y3NP*TANLC-Y6NP*TANUP+Y5PP*TAN0L>*( Y4PP*(-CR/2.~Y3NP*TANLC) 

2 + ( Y 4PP+Y5PP )/2.*(Y4PP#TANA- {— CR/2 Y3NP*TANLC ) ) — ( Y4PP**2 

3 4-Y4PP*Y5PP+ Y5PP**2) *TANA/3. ) 

IF I Y4PP.LT. Y5PP) XC BL£M= 

l ( CR+Y3NP*TANLC-Y6NP*TAN0P+Y4PP*TANGL)*( Y5PP*(-CR/2.-Y3NP*TANLC 
2 + Y4 PP# ( TANA-TANOE) > +( Y5 PP+Y4PP ) /2 .* ( Y5PP*TANDE— (-CR/2.-Y3NP+TANLC 
3+Y4PP*{ TANA-TANDE ) ) )— 1 . / 3.4= ( Y5PP**2+Y5P P*Y4PP+Y4PP**2 ) *TANDE ) 
CR=CR— CMDEX X 

S=2.*(-Y3NP**2*TANLC/2. + Y6NP**2*TAN0P/2.+YMIN*t CR+Y3NP*TANLC 
1— YfaNP*T ANGP+YMI N*TANOL/2 • ) + SADD + 

2 ( YMAX-Y MIN) /2.* (CR+ Y3NP*TANLC— Y6NP*TANQP+YMIN*T ANOL ) ) 

CR=CR+CHDEX X 


CBAR ,XC BAR, AND YCBAR ARE NOT VALID FOR A DISCONTINUOUS TRAILING EDGE 


YC8AR= 2 . /S * ( CR/ 2 .# YMIN **2+Y6NP **3*T ANOP/6.- Y3NP**3*TANLC/6 . 
1 +YM I N** 2*( YM I N/3.*T ANOL + Y3NP/2. *T ANLC— Y6NP/2»*T ANOP ) 

2+( YMAX+2.*YMIN) *( YMAX-YMIN) /6 . * ( CR+Y3NP *TANLC-Y6NP*TAN0P 
3+YMIN*TANOL ) ) 

C AV=S/ { 2 . *Y MAX ) 

IF( B IRA P- 82 RAP) 303,304,304 
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303 C8AR=2. /S*( CR**2*Y3NP+CR*Y3NP**2*TANPC*Y3NP**3*TANPC**2/3. + 

1(CR+Y3NP*TANLC) **2*( Y6NP-Y3 NP ) - ( Y6NP**2-Y3NP**2 ) * (CR+Y3NP*T ANLC ) * 
2TANLP-M Y6NP**3-Y3NP**3)*TANLP**2/3. + 

3 (CR+Y3NP*TANLC- Y6NP*TAN0P )**2*( YMIN— Y6NP ) +( YMIN**2— Y6NP**2I *T ANOL* 
4ICR+Y3NP*TANLC-Y6NP*TAN0P)+{YMIN**3-Y6NP**3)*TAN0L**2/3.+ 

5 { YMAX— YMIN) /3.* ICR+ Y3NP*T ANLC—Y6NP*TAN0P+YMIN*T ANOL )**2) 
XCBAR=2./S* (-YMIN/2.*CR**24-CR/2.*{ 1 . 5*Y3NP**2*T ANLC-Y6NP**2*T ANOP 
1/2.+YMI N**2 *< 3. *T ANA— TANE ) /2,-3.G*YMI N*Y3NP*T ANLC+YM I N* Y6NP * 
2TANOP)+Y3NP**3*TANLC*(2.*TANA-4.*TANC+TANP)/6.+Y6NP**3*TANA*TANOP/ 
36.+YMIN**3*TANA*TANQL/3.~Y3NP**2*YMIN*TANLC**2+lYMIN*Y3NP*Y6NP 

4— Y3NP/2 , *Y6NP J **2 ) *T ANOP*T ANLC+Y3NP/ 2 . *YMIN**2*T ANLC *( 2 .*TANA— T ANE ) 

5— Y6NP/2 . *YM IN**2*TANA*T ANOP + 

6XCBLEM ) 

GO TO 301 

304 CBAR=2./S*( CR**2*Y6 NP+CR*Y6NP**2*TANPC+ Y6NP **3* TANPC**2/3. +(CR 
1-Y6NP*TAN0P )**2*{ Y3 NP-Y6NP ) + { CR-Y6NP*TAN0P )* ( Y3NP**2-Y6NP**2 ) 

2*T ANOC+ ( Y3N P**3— Y6N P**3 ) *TANUC**2/3 . + 

3 { CR+Y3NP*TANLC— Y6NP*TAN0P )**2*(YMIN-Y3NP) + { CR+Y3NP*TANLC-Y6NP* 
4TAN0P)*TANDL*{YMIN**2-Y3NP**2) + ( YMI N** 3-Y3NP ** 3 > *TAN0L**2/ 3. + 

5 { YMAX-YMIN) /3.* (CR+ Y3NP*T ANLC-Y6N P*TANOP+YM I N*T ANOL ) **2 ) 
XCBAR=2./S*(-Y6NP/2.*CR**2+Y6NP**2/2.*( CR*TANC— CR/2 .*TANPC) +Y6NP** 
1 3/3 • *TANC*T ANPC — { Y3NP-Y6NP ) *CR/2 . * { CR-Y6NP*TAN0P )+( Y3NP**2- 

2Y6NP**2 ) /2. *( -CR*TANOC/ 2 . +T ANC*J CR— Y6NP *T ANOP ) ) + 1 Y3NP**3-Y6NP**3 i / 
33.*(TANC*TAN0C) + C YM I N-Y3NP ) * ( (CR+ Y3NP*TANLC-Y6NP*TAN0P) *!~CR/2 

4.-Y3NP*TANLC) ) + ( YMIN**2-Y3NP**2 )/2. *1 TANA*{ CR+Y 3NP*TANLC-Y6NP*T ANO 
5P)+TANOL*l-CR/2.— Y3NP*TANLC ).)■+( YMIN**3-Y3NP**3) /3.*TANA*TAN0L+ 
6XC8LEM ) 

301 ARN=4. * YMAX**2/ S 
AR6=ARN*BETA 
CR=CR/3 ET A 


DO 2 1=1 » NS MAX 
2 W(I)= 1 . 

DO 1301 I N= 1 t 10 
1301 BS ( IN 1 1 ) =0. 


READ IN THE 
THEY ARE 


VALUES OF THE COEFFICIENTS OF THE LOADING FUNCTIONS 
THE QP(JZ,l) 
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READ ( 5 * 511) ( QP( JZ» 1 ) » JZ=1,JKMAX) 


CHIB=ATAN(T ANC/ BETA ) 

ALAMB=ATAN( T ANA /BET A) 

PSI B = AT AN( T ANP/BETA ) 

OHE G8= AT AN( T ANE/BET A I 
TANC8=TAN{CHI8) 

TANLB^T AN( ALAMB) 

TANPB=TAN(PSIB1 
TANOB=T ANIOMEGB ) 

TANLCB=TANL B-TANCB 
TANLP8=T ANL B-TANPB 
TANPCB=TANP B-TANCB 
T ANGCB= TANOB— TANCB 
rANOLB=TANOB— TANLB 
TANOP8-TANO B-TANPB 
ETA ( I SS ST ) =0.0 
Y ( I SSST ) =0.0 
DO 7 NP=NNI I , NM AX 
ANP=NP 

ETA(NP) -SINI { ANP— SSSTA) *PI /(2.0*SSSTA) ) 

Y ( N P ) -E T A ( N P ) 

CHDSUB=0. 

I F ( E T At N P ) . GT « t Y6 NP— 0 .10) . AND. ETA (NP) .L E. Y6NP. AND.CHDEXX .NE .0. ) 
1CHDSUB=CHDE XX*( l.-t Y6NP-ETA (NP) 1/.10) /BETA 
IF( ETA(NP) . GT.Y6NP. AND.CHDEXX.NE.O. ) CHDSUB=CHDEXX/ BET A 
I Ft Y3NP.GE.Y6NP) GO TO 307 
IF( ETA( NP) . GE .Y3NP ) GO TO 309 
C ( N P ) — t CR+ETA ( NP ) *T ANPC B-CHOSUB ) / 2 . 0 
D(NP)=ETA{NP)*l TANP B+ TANCB ) /2. 0— CHD SUB/ 2 • 0 
DI FF3=F T A( NP ) — 0 . 

312 IF { ITTU.NE.2) GO TO 314 
I USX=NP— 1 
I UST=NP— 2 

C( IUSX)=tlO.O*C< I USX ) +2 .0*C ( I.UST ) ) / 12 . 

D ( I U S X ) — ( 1 0 • 0 *D t I US X ) +2 . 0*0 ( IUST 11/12- 
I TT U= 1 

314 IF ( 0 IFF 3.LT .0. . OR .0 IFF3.GT. .01) GO TO 323 
I TTU-2 

323 IF(NP-NNIi) 7,6,7 

309 IF(ETA(NP).GE.Y6NP) GO TO 327 

C ( N P ) = ( CR+Y 3NP*TANLC8— FT A (NP ) *T ANLP 8— CHDSUB ) /2 . 0 


78 



D<NP)=l-Y3NP*TANLCB+ETAlNP)*lTANPB+TANLB)-CHDSUB)/2.0 

DIFF3=ETA(NP)-Y3NP 

GO TO 312 

327 IFIETAI NPl.GE.YMIN) GO 10 331 

C ( NP) = ( CR.+Y 3NP*T ANLCB— Y6NP*TAN0PB+ETA IN P )#T ANOL B— CHOSUB ) /2» 0 
Dl NP )=( — Y3NP*TANLCB— Y6NP#TAN0PB+ETA (NP)*ITANQB+TANLB)— CHDSU8)/2.0 
DIFF3=ET A(NP)-Y6NP 
GO TO 312 

33i IFIYMIN.EQ.Y5PP) GO TO 340 

CATY4=CR+Y3NP*TANLCB-Y6NP*TAN0PB+Y4PP*T ANOLB— CHDSUB 
DAT Y4=X4PP/ BETA+CATY4/2 . 

C ( NP ) = . 5* ( CATY4— CATY4* (F.TA(NP)— YMIN)/(Y MAX— YMI N ) ) 
D(NP)=DATY4+(X5PP/BET A-DATY4 ) *IETA(NP)-YMIN)/(Y MAX— YMIN) 
DIFF3=ETAINP)-YMIN 
GO TO 312 

340 C(NP)=(Y4PP— ETA(NP) ) / 1 2 .0*1 Y4PP-Y5P P ) ) *1 CR+Y3N P#T ANLCB— Y6NP 
1*TAN0PB+Y5PP*TAN0LB-CH0SU3) 

DAT Y5=<-Y3NP*TANLC8-Y6NP* TAN0PB+Y5P P* ( T ANLB+TANQB ) -CHDSUB 1/2.0 
D(NP)-=DATY5 + (X4PP/BETA-DATY5)*( ETAl NP )- YMIN ) / l Y MAX-YM I N ) 
DIFF3=ETAINP)-YMIN 
GO TO 312 

307 IF C ETAl NP). GE.Y6NP) GO TO 311 

C(NP) = ( CR+ETA(NP)#T ANPC B— CHOSUB ) / 2. 0 
D(NP)=ETA (NP) *1 TANPB+TANCB) /2.0-CHDSUB/2.0 
DIFF6=ETA(NP)-0. 

313 IF C ITTU.NE.2) GO TO 316 
I USX=NP— 1 
I UST=NP— 2 

Cl IUSX) = (10.0*C( I USX ) +2.0*C ( IUST ) >/ 12 » 

D(IUSX)=(10.0*Dl IUSX)+2.0*D(IUST))/12. 

ITTU=1 

316 IFIDIFF6.LT.O..OR.DIFF6.GT. .01) GO TO 325 
ITTU=2 

325 IF(NP-NNII) 7,6,7 

311 IF ( ETAl NP ) . GE.Y3NP ) GO TO 329 

C(NP) = ( CR— Y 6NP*TAN0PB+ET A (NP ) *T ANOC B-CHDSUB ) /2. 0 
D(NP)=(-Y6NP*TAN0PB+ETA(NP) * 1 TANOB+TANC B) -CHDSUB) /2.0 
DIFF6=ETA(NP)-Y6NP 
GO TO 313 

329 IF( ETAl NPJ.GE. YMIN) GO TO 333 

C ( NP ) = ( CR+Y 3NP*T ANLCB— Y6NP* TANOPB+ETAlNP)*(T ANOLB )— CHDSUB ) / 2. 0 
D(NP)=( — Y3NP*TANLCB-Y6NP*TAN0PB+ETA(NP) *(TAN08+TANL B)-CHDSU8 1/2.0 
DIF F6=ET AINP)— Y3NP 
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n 

GO TO 313 

333 IF(YMIN.EQ.Y5PP) GO TO 346 

CAT Y4=CR+Y3NP*T ANLCB— Y6NP*TAN0PB+Y4PP*T ANQLB— CHDSUB 
0ATY4=X4PP/BETA+CATY4/2. 

C( NP )=. 5* (CAT Y4— CAT Y4*(ETA(NP)— YMIN)/ ( YMAX— YMIN ) ) 

D( NP )=0 ATY4 + ( X5 PP/BETA-OATY 4 )*( ET A( NP )— YMlN) /( YMAX— YMIN ) 

D I F F6=F T A ( N P ) — Y M I N 
GO TO 313 

346 C( N P ) = ( Y4PP— ETA (HP) ) / ( 2 .0*1 Y4PP— Y5PP) ) * ( CR+Y3NP*TANLCB-Y6NP 

1*TAN0P8+Y5PP*TANQLB-CHDSUB> 1 

0ATY5= (-Y3NP*TANLCB-Y6NP*TANQPB+Y5PP*(T ANLB + TANOB )— CHDSUB ) /2.0 
D ( NP ) =D AT Y5+ ( X4PP/B ETA-DATY5 ) *1 ET A( NPJ-YMI N) /( YMAX— YMIN ) 

D I F F6=E TA ( N P ) -YM IN 
GO TO 313 

6 C ( I SSST )=(5 .#CR+2.*C( NN1 I ) ) /12. 

D ( I S S ST ) =0 { NN 1 1 ) / 6 • 

7 CONTINUE 

DO 187 N-NN 1 1 tNMAX 

KSU=NMAX+1-N 

ETA (K SU ) =— E TAIN) 

Y ( K SU ) = — Y ( N ) 

C (KSU)=C { N) 

D(KSU)=DIN) 

187 CONTINUE 

CHI I=RAD*CHI 

ALAM=RAD*ALAMD 

ALADL=R AD*AL AME 

OHEG=R AD* OMEGA 

EOMG=RAD*EOMEG 

PS I 8B=R AD*P S I 

CR=CR*BETA 

CTB= ( X5-X4) /YMAX1 

TAPER=CT6/CR 

CFB=Y3NP*TANC 


WRITE INPUT DATA 


WRI 

WRI 

WRI 

WRI 


TE (6,501) 

TE(6,509) CASE, SYM, MACH 
TE ( 6 , 5 14 ) 

TE (6,527) C STA, SSSTA 
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WRITE16 

,510) 


WR I TE ( 6 

,526) 


WRI TE ( 6 

,555) 


WRI TE (6 

,532) 


WRI TE ( 6 

♦ 535) 

X2,Y2 

WRI TE ( 6 

,536) 

X3,Y3 

WRITE (6 

,537) 

X4,Y4 

WRITEI6 

,545) 

X5,Y5 

WRI TE ( 6 

,550) 

X6,Y6 

WRI TE ( 6 

,549) 

X6 A » Y6 A 

WRITE (6 

,556) 

X7,Y7 

WRI TE 16 

,510) 


WRI TE ( 6 

,558) 

YMAX1 

WRI TE ( 6 

,510) 


WRI TE (6 

,553 ) 


WRITE (6 

,555) 


WRITE ( 6 

,529) 


WRITE (6 

,503) 

X2PP,Y2PP 

WRI TE(6 

,504) 

X3NP, Y3NP 

WRITE 16* 

,505) 

X4PP,Y4PP 

WRI TE ( 6 

,506 ) 

X5PP,Y5PP 

WRI TE ( 6 

, 507 ) 

X6NP,Y6NP 

WRI TE (6 

,540) 

X6AP,Y 6AP 

WRI TE (6 

,508) 

X7PP, Y7PP 

WRI TE(6 

,510 ) 


WRI TE { 6 

,551) 


WRI TE (6 

,522 ) 

ARN, S, CAVfCBAR ,XCBAR,YC8AR 

WRI TE 16 

,548 ) 

CH 1 1 , ALAM, AL ADL , PS IBB, OM EG , EOMG , DELT , S I GM , RHH 

WRITE (6 

,524) 

CR , CTB , GFB , TAPER » BIRAP , B2RAP , XP , YP , Z P , CHDEXX 

WRI TE (6 

,502) 


MAIN PROGRAM 



00 79 K S= 1 * JMAX 
PKS=KS 

PHI (KS)=(2* G*PKS*PI )/(2*0*CSTA+l.0) 
XI KS ) =( I .—COS ( PHI (KS) ) ) /2. 

DO 79 NU=I SSST f NMAX 
CAN S 1=0 .0 

JK= ( KS-2 } SSST+NU+ 1 
ANU =NU 
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DO 14 N=l,NMAX 
AN=N 

VEIN) =COSl ( ( AN-SSST A)*PI)/I2,0*SSSTA) ) 

NNUD=IASS(N— NU) 

IF< NNUD.NE.O) GO TO 9 

BIN ) - (2 . O^SSSTA ) / (4.0*C0S { ( I ANU-SSSTA)*PI ) / ( 2.0*SSSTA1 I ) 

GO TO 14 

9 I F { MOD ( NNUD , 2 ) * EQ .0 ) GO TO 12 

8(N)=VE (N)/ ( ( 2. 0*SSSTA) *( ETAI N)— YINU) )**2) 

GO TO 14 
12 8 ( N ) = 0 . 0 
14 CONTINUE 

00 79 J = 1 , J MAX 
AJ = J 

DO 30 N=1,NMAX 

AK=0« 0 

AN=N 

IF(N.NE.NU) GO TO 16 
IF( J.EQ.l) GO TO 18 
IFIJ-2) 20,19,20 

18 AK=2.0*PHI( KS1+2.G* SINIPHIIKS)) 

GO TO 21 

19 AK= PH II KS ) - • 5 * S I N( 2 ,0*PH I ( KS ) ) 

GO TO 21 

20 GA= ( S I N( I A J-2 . 0) *{ PHI I KS ) ) ) )/ l AJ— 2.0 ) 

AK=GA- ( S IN ( ( AJ ) * (P HI( KS ) ) ) )/AJ 

21 PARTL=B (N) *AK 
A=0.0 

00 25 MUP—l , NMAX 
NUPNU=I ABS( NUP-NU) 

IF C NUPNU.EQ.O) GO TO 25 

IF I MODI NUPNU, 2 ) » EQ. 0) GO TO 25 

SSND=ABS(Y< NUJ-ETAI NUP) ) 

I F { SSND.EQ.O. ) GO TO 25 

ANUP=NUP 

UURR=ANUP— S SST A 

A= A + ( (COS( l UURR*PI ) / (2.0*SSSTA) ) )*# 2 ) *ALOGt SSND ) 

25 CONTINUE 

IF(J.NE.l) GO TO 28 

DF( 1)=--1.Q/(2.0*(SIN< IPHI(KS) )/2.0) )* I S INI ( PHI I KS ) > /2.0 ) > > 

GO TO 29 

28 OF I J )■=( A J— 1 . 0 ) * (COS ( I AJ— 1 .0 ) * I PH I IKS ) ) ) ) 

29 VL = 1,0/1 ((C(N>)**2 ) *2 ,0*SSSTA*V£( N)*S INI PHI (KS) ) ) 


82 



no noon 


AL! J, N) =VL*DF! J )*< C ,25*SSSTA)*( U Q~2, 0* t VE ! N > ** 2 > -ALGG ! 4. 0 ) )-A»+PA 
1RTL 

GO TO 30 


CHORDAL INTEGRATION SUBROUTINE 

SOLVES FOR THE CHORDAL INFLUENCE FUNCTION VALUES 


16 XSU8=-C(NU)*COS(PHI {KSH+D! NU) 

YSUB=Y ( NU) 

ETA$UB=ETA! N) 

CSUB=C i m 
DSU8=D( N) 

GO TO ! 351, 362, 353, 354,356,357,358*359* L070, 107 I) ,J 
1071 CALL GAUSS! Q., PI ,3 f SUMIO, FOFTIO ) 

AK= SUMl 0 
GO TO 355 

1070 CALL GAUSSI0.,PI,3,$UM9,F0FT9) 

AK= SUM9 
GO TO 355 

359 CALL GAUSStO* ,PI , 3 , SUM8 ,F0FT8 ) 

AK= SUMS 
GO TO 355 

358 CALL GAUSS! 0. ,P 1,3, SUM7 f F0FT7) 

AK= SUM7 
GO TO 355 

357 CALL GAUSS! G* ,PI ,2, SUM6,F0FT6J 
AK=SUM6 
GO TO 355 

356 CALL GAUSS! 0. ,PI ,2, SUMS ,P0FT5) 

AK= SUMS 
GO TO 355 

354 CALL GAUSS! 0. , PI ,2, SUM4,FQFT4) 

AK = SUM4 
GO TO 355 

353 CALL GA USS ! 0 . , P l , 2 , SU M3 f F OF T3 ) 

AR=SUM3 
GO TO 355 

352 CALL GAUSS! 0* y.PI.,.2, SUM2 f FOfTZ) 

AK= SUM2 
GO TO 355 

351 CALL GAUSS! 0* t P l « 2, SUMl , FOF T L ) 
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c 

c 


AK-SUM1 


355 AL( Jt N) “~B( N)*AK 
30 CONTINUE 

00 79 NP=IS SSTt NHAX 

1 = ( J-2)*ISSST+NP+1 

IF( NP.EQ.ISSST) GO TO 73 
NR=NMAX + l-NP 
I F ( ISYM.NE* l) GO TO 77 
SUHL « JK , I ) = AU J t NP ) + AL ( J , NR ) 

GO TO 78 

77 SUM .LI JK » I )~ AL( J , NP) - AL < J , NR ) 

GO TO 78 

73 I F ( ISYM.NE. L) GO TO 75 
SUHL ( JK f I )- AL(J f NP) 

GO TO 78 

75 SUMLCJK, I )= 0*00 000000 

78 CAN ST=C ANST +SUML < JK , I ) *QP ( I , 1 ) 
I F ( I .NE.JKMAX) GO TO 79 
CONST { jK f I.) -CAN ST 

YUT ( KS , .'NU) = -CON S T( JK, L) /4* 


IF(KS.NE.JMAX) GO TO 79 
CSP = 2 • *C ( NU )*B ETA 
IFC NU* EQ* ISSST) CSP=CR 
WRI T E ( 6 , 5 18 ) NU, ETA{ NU ) ,CS P 
WRI TE ( 6.9 510) 

WRITE (6 ,510 ) 

WRI TE (6.9.528) 

WRITE (6,521 ) 

WR I T E C 6 , 5 1 9 ) ( YU T ( KS UT , NU ) »K SUT= 1 , JMAX) 

WRITE (6,513) 

WRI TE ( 6 , 5 L6 ) < X ( KS ) KS= 1 , J MAX ) 

CALL LSQPOL ( X , YUT ( 1 ,NU) ,W ,RES ID, NS, SUMS T , LS , AST , BS, MS, CUT , NSMAX, MS 
1 MAX ) 

WRI TE (6 ,5 10 ) 

WRI TE ( 6 ,559) 

WRI TE ( 6 ,520 ) 

WR I T E ( 6 , 5 1 9 ) { 8 S ( J MS , 1 ) , J MS = 1 , 1 0 ) 

WRI TE (6 ,510 ) 

WRITE (6, 515) 
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u o 


WRITEI6 ,517 ) 

DO 700 I P T= 1 , i 1 
AI=IPT 

XC= (AI— l.)/10» 

CALL GAUSSt 1 . 0, XC » 1 , ZC » WU ) 
WRI TE (6 ,560 ) XC,ZC 
700 CONTINUE 


WRI TE ( 6 , 512 ) 
79 CONTINUE 
GO TO 1 
3 STOP 
END 
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FUNCTION F0FT1I THETA) 

COMMON XSUB »YSUB » ET ASUB ,C SUB, DSUB ,B S ( 10 , 1 ) , X ( 10) »XG,NS 
XI=-CSU8*C0S( THETA) +DSU8 

BK= 1 . + ( XSUB-XI I /SURT( (XSUB-XI )**2+( YSUB-ETASUB) **2) 
F0FT1=BK*(1 .0+C0S( THETA ) ) 

RETURN 

END 


FUNCTION F0FT2( THETA) 

COMMON XSUB,YSUB,ETASUB»CSUB,DSUB,BSl 10,1), XI 10),XC,NS 
XI--CSU3*C0S( THETA) +DSUB 

BK=l.+( XSUB-XI ) /SQRT( (XSUB-XI }**2+( YSUB-ETASUB ) **2 ) 

FOFT2=8K*Sl N( THETA) **2 

RETURN 

END 


FUNCTION F0FT3I THETA) 

COMMON X SUB, YSUB, ET ASUB, CSUB, DSUB, BS( 10,1), XI 10),XC,NS 
X I =— C SU B*CO S I TH E T A ) +DSUB 

BK=l.+( XSUB-XI) /SORT ( (XSUB-XI )**2+( YSUB-ETASUB ) #*2) 

EOF T3=BK*SI N( THET A) *S IN ( 2 .0 *THETA ) 

RETURN 

END 


FUNCTION F0FT4( THET A) 

COMMON XSUB, YSUB, ETASUB,CSUB,DSUB,BS( 10 , 1 ) , X ( 10 ) ,XC ,NS 
X I ~~C SU B*CO S ( THE T A) +DSU B 

BK=l.+( XSUB-XI) /SQRT( (XSU8-X I ) **2+( YSUB-ETASUB ) **2 ) 
F0FT4=BK*S l N( THETA) *SIN( 3.0*THETA) 

RETURN 

END 


FUNCTION F0FT5( THETA) 

COMMON X SUB, Y SUB , ETA SUB ,CSUB,DSUB,BS( 10 , 1 ) , X( 10 ) ,XC ,NS 
XI =— CSU8*C0 S (THET A) +OSUB 

BK= 1 • + ( XSUB-XI) /SQRT ( (XSUB-XI) **2+( YSUB-ETASUB ) **2) 
F0FT5=BK*SI N( THETA) *SIN( 4.0*THETA) 

RETURN 

END 



FUNCTION F0FT6( THETA) 

COMMON XSU8 » YSUB »ET ASUB »CSUB» DSUB »B SI 10, 1),X( 10),XC,NS 
XI=— CSUB*COS( THETA) +DSUB 

BK— l.+I XSUB— XI ) /SORT l (XSUB— XI )**2+( YSUB-ETASUB ) **2) 
FOFT6=BK*SIN(THETA)*SIN(5.0*THETA) 

RETURN 

END 


FUNCTION F0FT7I THETA) 

COMMON XSUB, Y SUB, ET ASUB ,CSUB, DSUB ,B S ( 10 , 1 ) , X ( 10 ) , XC ,NS 
XI=-CSUB*COS( THETA) +OSUB 

BK=l.+( XSUB-XI T/SQRTI (XSUB— XI )**2+( YSUB-ETASUB) **2) 

FOF T7=BK*SI N( THETA) *SIN (6,0*THETA) 

RETURN 

END 


FUNCTION F0FT8( THETA) 

COMMON XSUB, YSUB, ETASUB,CSUB, DSUB, 6S( 10, 1 ) , X (10 ) , XC ,NS 
XI =— CSU8*C0 S ( THETA) +DSUB 

BK- 1 . + ( XSUB-XI )/SQRT( (XSUB-XI )**2+( YSUB-ETASUB ) ** 2 ) 
F0FT8=BK*SI N( THETA) *SIN(7.0*THETA) 

RETURN 

END 


FUNCTION F0FT9( THETA) 

COMMON XSUB, YSUB »ETASU8,C SUB, DSUB »BS( 10,1 )*X(10)*XC, NS 
XI=-CSUB*COS( THETA) +DSUB 

BK=i.+( XSUB-XI ) /SQRT( (XSUB-XI )**2+( YSUB-ETASUB ) **2) 

FOF T9=BK*SI N( THETA) *SIN( 8.0* THETA) 

RETURN 

END 


FUNCTION FOFTIO (THETA) 

COMMON XSUB , YSUB ,ET ASUB ,CSUB, DSUB ,BS( 10 , 1 ) ,X( 10 ) »XC »NS 
XI--CSUB*C0S( THETA) +OSUB 

BK=l.+( XSUB-XI ) /SORT ( (XSUB-XI )**2+( YSUB-ETASUB ) **2) 

FOF T10=BK*S INC THETA ) *SI N ( 9. 0*THETA) 

RETURN 

END 




FUNCTION WU(XC) 

COMMON XSUB,YSUB,£TASJ8,CSUB,DSUB,BS1 10,1) ,X( 10) ,QQ,NS 
WU1=BS( i,l)+BS<2»l)*X{ 1 ) +8SI3 ,1 )*XC1 )**2+BSf 4, 1 )*XI 1 )**3+BS<5,i 
1)*X(1 )**4+BS(6,l)*X<l )**5+BSI7, i)*X| 1 ) **6+BSI 8, 1)*X| 1 )**7+BSI9 
2 » 1 ) * X ( 1 I ** 8+B S 1 10 , 1 ) * X ( 1 )**9 

WUN=BS{ 1, 1 ) +8 SC 2,1) *X(NS) +BSI 3 ,1) *XINS )**2>BS(4» 1 )*XINS)**3»BS C5# 1 

1 ) *X{ NS )**4+BS(6 , 1 )*X{NS1**5+BS ( 7 , l )*XI NS I**6+BS1 8, 1I*X(NS) ♦♦7+BSC 9 

2.1) *X(NS)**8+BS< 10, 1)*X(NS)**9 

DWU10X=BS ( 2, 1 ). +2 ,*BS(3 » 1 ) *X.C 1 )*3.*BS( 4, 1 >*X< l >**2+4.*BS< 5,1)*X(1 
1 )**3+5.*BS(6, 1)*X( 1 M'*4+6.*BSI7,i)*XU )**5+7.*BS 18, 1) *X ( 1 )**6 + 
28»*BSC9, 1)*X( 1 )**7+9.*BS< 10,1)*X< 1 )**8 
DWUNDX=BSI2,l)+2.*BS(3,ll*X(NS)+3.*BS(4,i)*X(NS)**2+4.*BSI 5, i)*X(N 
1S)**3+5.*BS(6» 1 )*X ( NS) **4*6. *8 SC 7, 1) *X C NS) **5*7, *8SC8,i)*X(NS)**6+ 
28»*8S(9»1)*X(NS) **7+9 ,*BSI 10, 1 ) *X(NS ) **8 
IFIXC.LE.Xt 11 .AND.XC.GE.O.) WU®WUI+DWU1DX*(XC-X< 11.) 

IF (XC.GT .X ( 1 ) .AND.XC.LT «X(NS ) ) 

1WU = BS<1»1M-BS{2,1)*XC +BS(3,1)*XC **2+BS(4, 1)*XC **3 + BS(5,l 

2 ) *XC **4 +B S ( 6 , 1 1 *XC **5+BS I 7* 11 *XC **6+BS( 8, 1)*XC **7+BSI9 

3.1) * XC **8+BS{ 10,1)*XC **9 

IF C XC.GE.X(NS) . AND. XC.LE. 1,1 WU=WUN+DWUNDX*1XC-X<NS) 1 

RETURN 

END 
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o o o non n n n o o o 


SUBROUTINE MAT I NV I A ,N, 8 , M, DETERM, IP IVDT , INDEX, NMAX, ISCALE) 

MATRIX INVERSION WITH ACCOMPANYING SOLUTION OF LINEAR EQUATIONS 

DIMENSION IPIVOTCN) , AC NMAX, N) , B ( NMAX, M) , INDEX (NMAX, 2) 
EQUIVALENCE C IROW,ORQW>, ( I COLUM , JCOL UM ) , CAMAX, T, SWAP) 

INITIALIZATION 

5 I SC ALE-0 

6 R1=10.0**100 

7 R2-1 • O/Rl 
10 Q£TERM=1.0 
15 DO 20 J= L »N 
20 IPI VO T i J ) =0 
30 DO 550 1=1 , N 

SEARCH FOR PIVOT ELEMENT 

40 AM A X= 0.0 
45 DO 105 J= 1 , N 

50 IF I IPI VOTC J)-l) 60, 105, 60 
60 DO 100 K=lfN 

70 IF CIPIVOnK)-!) 80, 100, 740 

80 IF ( A8S ( AMAX ) — A8S ( A ( J , K ) ) )85,100, 100 

85 I ROW= J 

90 ICOLUM=K 

95 AMA X= A ( J , K ) 

100 CONTINUE 

105 CONTINUE 

IF CAMAX) 110,106,110 

106 D£ TERM- 0.0 
I SC A L 6=0 
GO TO 740 

110 IPI VO T ( I COLUM) = I PIVOT C I COLUM ) +1 

INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL 

130 IF ( I ROW- I COLUM ) 140, 260, 140 
140 DETERM=-DET ERM 
150 DO 200 L=i , N 
160 $WAP= A( I ROW ,L ) 
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u o u o o 


170 
200 
205 
210 
2 20 
230 
250 
260 
270 
310 


1000 

1005 

1010 


1020 


1030 

1040 


1050 

1060 

1070 


1080 


1090 

2000 


2010 

320 



:: 

| 


A( I ROW , L ) =A ( ICOLUM, Li 
A( I COLUM, L) — SWAP 
IF ( M I 260, 260, 210 
DO 250 L=1 » M 
SWAP— B( I ROW » L ) 

B( I ROW, L )=B I ICOLUM, L ) 

B( ICOLUM, L)=SWAP 
INDEX ( I , 1 )= I ROW 
INDEXII ,2)=IC0LUM 
PIVOT=A( ICOLUM, ICOLUM) 

SCALE THE DETERMINANT 

P I VUT I = PI VO T 

IF C ABSI DETERM )-R 1)1030, 1010, 1010 
DETERM=DETERM/R 1 
I SC ALF= I SCALE+ 1 

IF (ABSI DETERM)— R1 ) 1060, 1020,1020 
DETER M=DETERM/R1 
1 SCALE= ISCALE+1 
GO TO 1060 

IF I ABSI DETERM )-R2) 1040, 1040, 1060 
DETERM=DETERM*R1 
I SC ALE = I SCALE— 1 

I F ( ABS( DETERM )-R2) 1050, 1050,1060 
DET ERM=DETE RM*R 1 
I SC AL E= I SCAL E— 1 

IF t ABSI PIVOT I )-Rl ) 1090, 1070,1070 
PIV0TI=PIVCTI/R1 
I SCAL E= ISCALE+1 

IF! ABSI PI VOTI )-Rl ) 320, 1080, 1080 
PIV0TI=PIV0TI/R1 
I SC ALE= I SCALE+ 1 
GO TO 320 

IF! ABSI PIVOT I )-R2) 2000, 2000,320 
P I VOTI=P l VOTI #R 1 
ISCALE= 1 SCALE— 1 

I F ( ABSI PI VOTI )-R2) 20 10, 20 10, 320 
PIVOT I=P I VOT I*R l 
I SC ALE= I SCALE-l 
OETERM=DETE RM*P I VOT I 

DIVIDE PIVOT ROW 6V PIVOT ELEMENT 
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c 

330 A( ICOLUM, ICOLUM) =1» 0 
340 DO 350 L=1,N 

350 A( ICOLUMfL) =AI ICOLUM»L) /PIVOT 
355 IFIM) 380, 380, 360 
360 DO 370 L=1,M 

370 B( ICOLUM»L)=BI ICOLUM»L) /PIVOT 
C 

C REDUCE NON-PIVOT ROWS 

C 

380 DO 550 Ll=l,N 

390 I F C LI— ICOLUM) 4 00, 550, 400 

400 T=AIL1, ICOLUM) 

420 A( L 1 , ICOLUM ) =0. 0 
430 DO 450 L=1,N 

450 A{L1,L)=A(L1,U-A(IC0LUM,L>*T 
455 IFIM) 550, 550, 460 
460 DO 500 L=1 * M 

500 B(LI,L}-8ILl,Ll— BlI COLUM»L) *T 
550 CONTINUE 
C 

C INTERCHANGE COLUMNS 

C 

600 DO 710 1=1, N 
610 L=N+l-I 

620 IF ( INDEX I L , 1 )- I NDE X ( L, 2 ) ) 630, 710, 630 
630 JR0W=INDEXIL,1) 

640 JCOLUM= I NOE X I L , 2 1 
650 DO 705 K= 1 , N 
660 SWAP=AI K , JROW) 

670 AIK, JROW)=A(K,JCOLUM) 

700 AIK, JCOLUM) =SWAP 
705 CONTINUE 
710 CONTINUE 
740 RETURN 
END 

CANE206 JEAN MIGNEAULT 

C LEAST SQUARE POLYNOMIAL FIT ANE206 
C 



SUB ROUTINE L SQP OL < X , Y , W , RES I D , N , SUM , L , A , B , M , C » N MAX , MMAX 1 

DIMENSION X ( NMAX) , Y { NMAX » LI , RES I DINMAX, LI # AIMMAX»MMAXI » 

2 B I MM AX * U . C ( N MAX , Ml » SUM ( L ) , W ( NM AX ) 

10 00 20 1=1, N 

20 C( I ,11=1.0 

30 DO 50 J=2,M 

40 DO 50 1 = 1, N 

50 CU,J1 = C< I, j-il*X(I 1 

.60 DO 100 1 = 1, M 

70 DO 100 J=1 » M 

80 A( I , J ) =0.0 

90 DO 100 K= 1 , N 

100 A( I ,J) = A( I, J1+C(K,I )*C(K, J)*W(K) 

105 DO 150 J=1 , L 
110 DO 150 1=1, M 
120 8( I ,J)=0.0 
130 DO 150 K=1 , N 

150 B ( I , J 1 = B ( I , J 1 +C t K , I ) #Y ( K , J ) *W ( K ) 

CALL MATINV I A, M,B, L, DETERM , RES ID, C,MMAX, ISCALE1 
180 DO 205 4=1, L 
185 SUM ( J ) =0. 0 
KK=M 

192 DO 195 K=1 , M 
C (K » 1 ) = B< KK , J ) 

195 KK=KK-1 
198 DO 205 1=1, N 

RES I D ( I , J 1 = P0LYE1 (X ( 1 1 , M,C) -YI 1,41 
205 SUM ( J ) = SUM( J 1 +RESIDI I , J 1 **2*W( 1 1 
210 RETURN 
END 
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FUNCTION POLYEl 1 X »M »C ) 
DATA BIG/Q3 77 177777 777/ 
DIMENSION CIM) 

I F ( M— 1 ) 10 » 1 1 » 12 
12 N=M— 1 

P0LYE1=C( 1 ) 

D020 1=1 *N 

20 POL YE 1= X* POLYE.l +C ( I +1 ) 
RETURN 

10 POL YE 1= BIG 
RETURN 

11 P0LYE1=C(1) 

RETURN 

END 
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SUBROUTINE GAUSS ( A , B, N , SUM , FOFX ) 

REFERENCE SCARBOROUGH NUM. MATH. ANAL. PAGE 147 
HOWEVER THIS SUBROUTINE INTEGRATES FROM ZERO TO ONE 

DIMENSION U(5),R(5) 

U(l) =.425562830 5091 84 
U(2) =.283302302935376 
U( 3 )*. 16029 5215 8504 88 
U(4)=. 06746 83 16655508 
U( 5) =.0130467357414 14 
R 1 1 1 — . 14776 2112357376 
R I 2)=. 134633359654998 
R <3 )=. 109543181 257991 
R (4 )=.0 7472 56745 75290 
R< 5 )=. 03333 5672 154344 
SUM-0 • 0 

IF(A.EQ.B) RETURN 
F INE=N 

DELTA=F INE/IB-A) 

DO 3 K= 1 * N 
X I = K— 1 

F INE=A+X I /DELTA 
DO 2 11= 1,5 
UU=UI II ) /DELTA+FINE 

2 SUM=R (I I ) *FOFXC UU ) + SUM 
DO 3 L=l,5 

UU=(l.0-U(L)) /DELTA +FINE 

3 SUM=R I L ) *FOFX I UU) +SUM 
SUM= SUM/DELTA 
RETURN 

END 
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Sample Output Listing 
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STATION NUMBER" 29 SPANWISE LOCATION" .64945 


o 

• 

H 


8 
o o 

O. 

o 

p-4 O 

• 

V. 

< o 

2- 

a 

a 

X 

• 

o 



O' 


o 




• 


o 






o 




& u 


O' o 




• v 


c o 




2 X 


• 




O 


o 




u 

o 






•X 






X 





00 

uu 

o 




.* 

o 

o 






o 




a o 

on 

00 O 




* 

a c 

< o 




2 X 

UJ 

• 




a 

2 

o 



QC 

o 

o 




< 


a 




tu 


o 

2 




OC 





o 

?*»• 

1-4 

o 



h- 

i 

on 

o 





< 

o 




a o 

UJ 

f- o 



2 

t s 

c£ 

< o 



O 

2 X 

U 

* 



oc 

O 

2 

o 



u. 

o 

*“* 




X 


UU 




o 


o 




a£ 

uu 


a: 

o 



» 

• 

UJ 

o 



on 

H 

2 

h 

0- U 

* s 

O 

a: 

o 

o 
>0 o 
< o 

UJ 

Cl 

< 

(MOOsOCNjOaoOOO'O 

43 ® ^ ^ n «t in <n h o 

cjcvjsoinsofooo^oooocsjo 

Ip* 

2 X 

2 

o 

X 

voo^h^O^cO'OO'^O 

o 

a 

O 

o 



on 

i44-h-aors->0't< , vJ— <ooo 

«j 


» 


QC 

o 

o 


U 


UJ 


QC 

m 

N* 

>«; 

o 

X 

OOOOOOOOOOO 

£ 


o 

< 

ooooooooooo 

2 

• 



o 

ooooaooooooo 

o 

o 

H- 

a o 

• "s, 

U. 

o 

in o 
<t o 

§ 

UJ 

'vOOOOOOOOOOo 

XOHMiri'tlft'ONooff'O 

h* 

< 

2 X 

O 

co 

• 

o 

o 

•* 

u 

1- 


X 




2 






UJ 




*«s. 


Up 




3 f 


U 





4 - o 

pp 

O' 




• o tn 

u. 

4 



oo 

v- oo 

LL 

hr 



lii 

a m o ^ 

UJ 

4 - in 



Ql 

t 4 "N. o 

o 

< o 



o 

2 • X O' 

u 

• 



*1 

a i • 


in 



on 

o 

*1 

(M 





< 

1 





NP 






X 





ro sO 

i 

m 




* O' o 

>• 

in 




h CO O 

_J 

m 




Cl K u o 

o 

cn r- 




• r- v* tn 

a 

< s0 




2 • x r- 


• 




O 1 • 






o 


in 




r\J pH 

rsl 

* O 03 

h- 

U- in pH 

4 

a cn u 

(M nj 

• >0 'N ^ 

4* 

2 • X 4- 

» 

CD p-4 • 

*H 

O 1 



pH m 


4 

• m 

00 

m 

u- in 

cr 

4" 

a oo 

U 'O 

— 1 O' 

• (\J 

^ pH 

< CvJ 

2 • 

X pH. 

• 

o -* 

* 

4 


107 



STATION NUMBER* 30 SPANWtSE LOCATION* .70711 CHORD* 



o 


I— 

a o 

* v 

2 X 

a 

u 


o 
o 
o o 
o 
< o 


cn 

on 

o 



© 


O 


• 


O 


t- 


© 


a u 


O' o 


• \ 


< © 


z X 


* 


8 


© 



O 








X 



co 

• 

OF 

o 

o 


K 


o 


flL U 

00 

oo o 


• S 

oc 

< © 


2 X 

UJ 

• 


O 

X 

o 

a: 

o 

O 


< 


a. 


UJ 




a: 


o 




2 


O 

h- 

*« 

o 

K* 

t 

on 

o 


H* 

< 

o 

K 

a. u 

UJ 

r*- o 

Z 

* v 

a: 

< o 

o 

Z X 

o 

• 

OC 

o 

z 

o 

u. 

o 

*-* 


X 


IL 


o 


a 


0C 




UL 

>0 

a: 

o 


• 

UJ 

o 

CO 

L- 

o 

o 

K 

a o 

oc 

s 0 o 

X 

• N. 

o 

< © 


2 X 


• 

o 

O 

2 

© 

a 

o 

M 






-j 




a 


o 


oc 


•V. 


fe 

in 

X 

c 

X 

• 


o 

o 

L- 


o 

u 

a. o 

u. 

in o 



o 

< © 

K 

2 X 


• 

< 

O 

co 

o 

*> 

o 

L- 


** * 


2 


O 


UJ 




M 


3S 


o 


*** 

•4- ^4- 

»— 1 

00 

* 

• f— in 

CL 

00 

CO 

H* O 00 

UL 

O' 

UJ 

o- m o O' 

UJ 

-4* <4- 

CL 

• '4' N ■■* 0 

o 

<c m 

O 

Z • X O' 

o 

• 

-J 

O 1 • 


o 


u 

-J 

IM 



< 

1 



M 




X 




a 



CO rH 

2 

r- 


• CO o 

> 

in 


1 - 00 o 

-J 

in 


a. co o o 

o 

on in 


• m 

a. 

< co 


2 • x r- 


• 


o 1 • 




u 


in 


nj O' 

© 

* *-* co 

in 

W . f»4 

m 

a >0 um 

CM CsJ 


< m 

z • x >r 

• 

O H • 

CO 

© 1 

fO 


I 


uj m in r*- ir> -4 f\j (\i o 

O. 

l/> N^MB^'O'tNHOOO 


o£ 

UJ 

CO 


< 

o 


§ 


o 


ooooooooooo 
ooooooooooo 
u ooooooooooo 
^ooooooooooo 
XO*-*<v4<n'tinOf-0QCT'O 

O -l 


CO 

O' 

» m co 

m 

K "4“ © 

r- 

Cl co O <0 

O' 

• S. ~4 

< o 

Z • X ~4 

• 

o • 

u 



108 



STATION NUMBER= 31 SPANWISE LOCATION^ •76041 CHGRD= .50094 
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0N= .99692 CHORD= .12380 
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H. SUPPLEMENTARY PROGRAMS (A1590 AND A1591) 


ASPECT RATIO PROGRAM (A1590) 


To find the aspect ratio and root chord of a planform that is of the general type 
shown below (and which does not contain a trailing- edge inboard chord- extension) requires 
a knowledge of most of the quantities shown in figure 2. 



The term ARRATO is set equal to zero unless a particular value of aspect ratio is 
sought. If set equal to a value other than zero, the program then fixes the leading-edge 
sweeps, leading- and trailing- edge breaks, overall taper ratio, and tip sweepback of the 
planform, and iterates on the trailing-edge sweep angles (which have been initially speci- 
fied to be lower than an assumed final value) until the desired value of aspect ratio is 
obtained. This procedure works best if both trailing-edge sweeps are the same since they 
are adjusted together. 
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The input data are to be read in a 6F12.5 format where the quantities CHI, LAMBDAJ 
PSI, OMEGA, Bl, and B2 are on the first card and TAPER, BELT, and ARRATO are on 
the second card. 

The sample cases included are 

| 

1st: The planform like that shown in figure 2 

2d: A sweptback planform with skewed tip that requires use of the iteration 

procedure 

3d: A double delta type planform with sweptforward trailing edge 

4th: An arrow-type planform with skewed tip to be used in the pivot determining 
program 
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PIVOT DETERMINING PROGRAM (A1591) 


If the high sweep position of a variable -sweep wing is known and its planform has 
no leading- or trailing-edge break or trailing-edge chord-extension, like the planform 
shown in figure 3, and its subsonic aerodynamic characteristics or mean camber surface 
are required at this or lower sweep angles, it is necessary to first know the wing in its 
outer panel streamwise tip position in order that they may be computed by using Langley 
computer program A0313 or A0457. There are, however, an infinite number of these 
outer panel streamwise tip positions that will result in the same high sweep wing given 
the proper pivot location. They may be thought of as having been generated by a combi- 
nation of two variables: (1) ratio of low to high sweep semispan, and (2) the fractional 
location of the pivot along a chord which is normal to the high-sweep leading edge (T). 

(See fig. 4.) 

This program determines the absolute location of the pivot point relative to a new 
coordinate system and the resulting low sweep streamwise tip planform for the conditions 
given in the preceding paragraph. (See sample listing.) In case a pivot is selected which 
would result in the trailing edge inboard and trailing edge outboard not intersecting on the 
right side of or at the plane of symmetry in the streamwise tip position, the program 
reduces the root chord by an amount equal to 

CR(1 - X) + 0.0001(tan xfs - tan Q,) + tan 0 - tan A - Bl(tan % - tan A) 

(symbols are defined in fig. 2) and calls this the amount of trailing-edge chord- extension 
that must be specified with the streamwise tip position in those programs. The num- 
ber 0.0001 is assigned to be the new trailing-edge break location, and B2 is the old value. 
Correspondingly, the taper ratio is changed from being the actual skewed tip length over 
the root chord in the high sweep position to the same tip length but over the new root 
chord. The x- location of the pivot, as shown on the output sheets, is relative to the orig- 
inal half root chord but scaled to the new streamwise tip semispan so that it can be used 
for input to program A0313 or A0457 directly even if the root chord has been reduced by 
a trailing-edge chord- extension. 

The required data to be used for input are associated with the high sweep wing 
and are punched according to a 6F12.5 format. The first card contains the aspect ratio, 
leading-edge sweep angle, trailing-edge sweep angle, taper ratio, tip skew angle, and 
root chord for a unit high sweep semispan. On the second card are the values for span 
increase and maximum number of pivot locations desired. The last card(s) contains the 
table of fractional pivot locations whose maximum number has been specified on the 
second card. 
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Sample cases are given in the output listings for the planform shown in figures 3 
and 4. The T = 0.75 case is used in program A0313 as the first sample case with 
A = 45° and M = 0.60. 




I I 1 1 1 

O 4 .8 L2 1.6 


Figure 3.- Arrow planform in original high sweep position. 
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ASPECT RATIO PROGRAM A1590 
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Program Listing 
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PROGRAM ASPRATO (INPUT , OUTPUT , TAPE5= INPUT , TAPE6=0UTPUT) 

REAL LAMBDA 

500 FORMAT ( 6F12.5 ) 

501 FORMATt 11F9.5) 

502 FORMATt 1H1 , 4X2HAR6X3HCH I4X6HL AMBD A6X3HPST 3X 5H0MEGA7X2HB17X 
1»2HB2»5X»5HTAPER»4X«5HDELTA»6X'» 2HCR »7X * 2HB3 ) 

WRITE (6 .502 1 

1 RE AD ( 5 » 500) CHI , L AM8DA , PS I , OMEGA, B1 , B 2 
I Ft EOF » 5) 21 » 20 

20 RE AD (5.500) TA PER .DELTA , ARRATO 
ITEST=1 
JTEST = 1 

CHII=CHI/57. 29578 
AL AMB-L AMBDA/ 5 7. 2 9578 
DELT=-DELT A/57. 29 578 
TANC = SIN( CHII )/ CO St CHI I 1 
T ANL=SI Nt ALAMB) /CCS t ALA MB) 

TAN9D=SIN( 1 .57G7963+DELT ) /COS t l. 5707963+DELT ) 

3 IFtPSI. GT. LAMBDA. AND. TAPER. LE. 1.00) GO TO 15 
IFt PS I. LE. LAMBDA. AND. TAPER.LE. 1.00) JTEST=1 
GO TO 16 

15 JTEST=JTEST+1 

IF ( JTEST. GT .2 1 GO TO 1 

16 PSII=PSI/57. 29578 
0MEG=CMEGA/57.29578 
TANP=SIN(PSII )/COStPSII ) 

TANO=SI NtOMEG) /COS(OMEG) 

CR= 1 . / 1 l.-T APER*COSt DELT ) *t l.+TANO/T AN9D) ) *( Bl* t TANC-T ANL) +B2*( 
1TAN0-TANP) +TANL-T ANO) 

B3=1.-CR*TAPER*C0S(DELT ) /TAN9D 

S =2. * (-81 **2* tTANL— TANC )/2.+B2**2*(TANO-TANP)/2.+B3*(CR+Bl*tTANL- 
1TANC)— B2*(TANQ-TANP)+B3*(TAN0— TANL )/2. ) + ( l.-B3)/2.*tCR+Bl* 

2 IT ANL-T ANC )— B2*t T ANO-T ANP) + B3*l TANO-T ANL) ) ) 

AR = 4./S 

IFtARRATO.EO.O.O ) GOTO 2 
IF tABSt ARRATO-AR) .LT.O. 00001) GO TO 2 
IFt ARRATO. GT.AR. AND. ITEST.EO.l) GO TO 4 
IF (ARRATO.LT . AR. AND. ITEST.EQ.l ) GO TO 5 
IFt ARRATO. GT . AR. AND. IT EST . EO. 2) GO TO 6 
I F (ARRATO.LT .AR.AND.ITEST.E0.2) GO TO 7 
IF(ARRAT0.GT.AR.AND.ITEST.EC.3) GO TO 8 
IFt ARR A TO. LT. AR.AND.IT EST .EQ.3) GO TO 9 
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IFIARRATO.GT.AR. AND. ITEST.EG.4) GO TO 10 
IF (ARRATO.LT. AR. AND. ITEST.EQ.4) GO TO 11 
I F ( ARRATO. GT.AR.AND.1T EST . E0.5 ) GO TO 12 
IF ( ARR ATO .LT . AR. AND. ITEST.EQ.5I GO TO 13 
IFtARRATO.GT. AR. AND. IT EST . EQ.6I GO TO 14 

4 0MEGA=0PEGA+10 * 

PSI=PSI+10. 

GO TO 3 

5 OMEGA=CMEG A-9. 

PSI=PSI-9. 

ITEST=2 

GO TO 3 

6 0MEGA=CMEGA+1 . 

PSI=PS!+1. 

GO TO 3 

7 OMEGA=CMEGA-0.9 
PSI=PSI-0.9 

I TEST =3 
GO TO 3 

8 0MEGA=0 MEGA+O • 1 
PSI=PSI+0.1 

GO TO 3 

9 OMEGA=CMEGA— 0. 09 
PS I=PS 1-0 .09 
ITEST=4 

GO TO 3 

10 0MEGA=GMEGA+0. 01 
PSI=PSI+0.01 

GO TO 3 

11 0 MEG A=C PEG A— 0.009 
PS I =P SI -0.009 

I TEST=5 
GO TO 3 

12 OMEGA=CMEGA+0.001 
PSI=PSI+0.001 

GO TO 3 

13 0MEGA=CPEGA-0. 0009 
PSI=PS 1—0.0009 

I TEST=6 
GO TO 3 

14 0MEGA=CMEGA+0 . 000 1 
PSI=PSI+0.0001 

GO TO 3 
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IHB3.LE.0.) C 7 O TO 1 

WRITE (6, 501 ) AR, CHI , LAMBDA , PS I , OMEGA, B 1 ♦ R2, TA PER, DELTA ,CR, B3 
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PIVOT DETERMINING PROGRAM A1591 
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Program Listing 
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PROGRAM PIVOT ( INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT) 

DIMENSION XCR (1 00 )»TANL 0(100) ,TANLLO( 100) , T ( 100} 

REAL LAMBDA, LAM 
500 FORMAT! 6F 12. 5 1 

550 FORMAT! 1H1 ,40X, 10HINPUT DATA////} 

551 FORMAT! 25X,50HARRQW WING WITH CROPPED TIPS (HIGH SWEEP POSITION}} 

5 52 FORMAT ( 3X »3HAR = *F9. 5 , 1 X »4HCHI= , F9 .5 , 1 X ,4HPS 1= , F9. 5, IX , 6HTAP ER*, F9 • 

15,1X,6HDELTA=,F9.5,1X,3HCR=,F9.5, 1X,18HSPAN INCRSE DESR0=,F9.5) 

553 FORMAT! 46X,17HWING TIP LOCATION) 

554 FORMAT! 34X,9HX LE TIP=, F9.5,5X, 18HY LE TIP* 1 . 00000/34X,9HX TE TI 
1P=,F9.5,5X,18HY TE TIP= 1.00000////) 

555 FORMAT! 34X,9HX LE TIP* ,F9. 5,5X, 9HY LE TIP=,F9. 5/34X.9HX TE TIP=,F9 
1.5,5X»9HY TE TIP=,F9.5//) 

556 FORMAT! 20X,70HX AND Y LOCATION AND SLOPE OF LOCI OF ALL PIVOT POIN 
ITS THAT WILL SWEEP) 

557 FORMAT! 20X,68HTHE HIGH SWEEP WING OUTER PANEL INTO THE SAME LOW SW 
IEEP AND SEMISPAN//) 

558 FORMAT (47X»8HX P I VOT=, F9.5/47X, 8HY PI V0T=,F9. 5/49X, 6HSL0PE=,F9.5// 

l/> 

559 FORMAT! 10X, 99H*T*,FRACT ION OF CHORD NORMAL TO THE LEADING EDGE OF 

1 ARROW WING WHERE THE PIVOT IS SPECIFIED TO ACT,/9X,100HAL0NG WITH 

2 THE CORRESPONDING X INTERCEPT AT ROOT AND THE SLOPE OF ALL POINTS 

3 WHICH HAVE THE SAME *J*//) 

560 F0RMAT(35X,20HT, FRACTION OF CHORD*, F9.5/31X ,24HX AT ROOT AFT OF OR 
II GIN* » F9. 5/31 X,24HSL0PE THROUGH X AT ROOT*, F9. 5/// ) 

561 FORMAT! 1H150X, 11HCUTPUT DATA) 

562 FORMAT!///// 45X,26HWING IN LOW SWEEP POSITION/20X,71HALL OF TH 
IESE DIMENSIONS ARE SCALED FOR A UNIT SEMISPAN IN THIS POSITION//// 
2) 

563 FORMAT! 2X »3HAR* »F9. 5,5 X »4HCHI = *F9 . 5»5X»7HLAMBDA=»F9.5»5X»4HPSI=»F9 
l.5»5X6H0MEGA*»F9. 5,5X1 1HY LE BREAK=,F9.5/2X11HY TE BREAK*, F 9. 5, 5X, 
26HTAPER*,F9.5, 5X , 3HCR* ,F9.5,5X,31HSEMISPAN OF WING IN HIGH SWEEP*, 
3F9.5/50X, 16HCH0RD EXTENS I0N=,F9.5 ) 

564 FORMAT! 5X99HALL DIMENSIONS REFERENCED TO COORDINATE ORIGIN AT TNE 
1HALF ROOT CHORD AND SCALED ON A UNIT SEMISPAN///) 

565 FORMAT! 15X83HPIV0T LOCATION WHICH OCCURS AT THE INTERSECTION OF TH 
IE LOCI OF ALL PIVOT POINTS AND/ 17X6 1HTHE LINE WHICH CONTAINS ALL 0 
2F THE *1* FRACTION NORMAL CHORDS/ 47X, 8HX PI V0T=,F9. 5/47X,8HY PIVOT 
3=,F9.5/////) / 

566 FORMAT!///) 

567 FORMAT! 44 X,21HIN LOW SWEEP POSITION) 

568 FORMAT! IH1) 
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o o no 


1 REAOC 5, 500) AR ,CH I, PSI , TAPER, DELT A,CR 
IF( EOF » 5 ) 5,3 
3 READ( 5,500)SI,TMAX 
ITMAX=TMAX 

READ! 5, 500) ( TIN ) »N=1 , l TMAX) 


X4= 1.00000 
X 5= X4 + T A P E R *C R 
WR ITE ( 6,5501 
WRITE (6 ,551} 

WRITE ( 6,564) 

WRITE (6,552) AR,CHI,PS I , TAPER, DEL TA,CR,SI 

L AMBD A=CH 1/57 . 295 78 

0MEGA=PSI/57. 29578 

DELT A=DELTA/57. 29 578 

LAM=( LAMBDA-DELTA 1*57.29578 

OME= ( OMEGA -DELT A )*57. 29578 

TANC=SIN( LAMBDA) /COSILAMBOA) 

TANL=S I N( LAMBDA-DELTA) /COS (LAMBDA-DELTA) 

TANP= SIM OMEGA) /CCS (OMEGA) 

T ANO=S I N( OMEGA-DELTA) /COS (OMEGA— DELTA) 

T AND=S I N( DELTA) /CCS (DELTA) 

TAND2=SIN( DELTA/2.) /COS (DELTA/2.) 

Y4PP= 1.00 

Y5PP=(— CR + Y4PP*(T ANC+1 • /TAND) )/ (T ANP+1. /TAND) 

X4PP=-CR/2.+Y4PP*TANC 

X5PP=CR/2.+Y5PP*TANP 

WR I TE ( 6 ,5 55 ) X4PP ,Y4PP , X5PP,Y5PP 

A H=X4PP-X4 

A 12=Y4PP— SI 

A2l=X5PP-X5 

A22=Y5PP— SI 

C 1=0. 5* ( Y4PP**2+X4PP**2— SI**2-X4**2 ) 

C2=0. 5* l Y5PP**2+X5PP**2— SI **2-X5**2) 

DETERM=A11*A22-A12*A21 

XP=( l./DETEPM)*( A22*C1-A12*C2) 

YP=(1./DETERM)*(-A21*C1+A11*C2) 

WR I TE ( 6 ,556) 

WRITE (6 ,557) 

WRITE ( 6,558) XP,YP,TAND2 



WRITE (6*567) 

WRITE! 6*555) X4,SI»X5,SI 
WRITE (6 »56II 
DO 4 1 = 1* ITMAX 

TANLO ( I )=S IN ( ! L AM BDA-OMEGA )*T (II) /COS ( ( LAMBDA-OMEGA )*T ( I ) I 
XCR! I I = (T ( I )*COS ( LAMBDA l+S!N(LAMBOAl*T ANLO! I) » *CR/ (COS (LAMBDA! 
H-SIN! LAMBDA)*T ANLO( III -CR/2. 

TANLLO( I) =SIN (LAMBDA- (L AMBCA-OMEG A)*T (III /COS ( LAMBDA- ( LAM8DA- 
10MEGA)*T( III 
WRITE (6 *562) 

WRITE! 6*559) 

KfPsYfR I I I /5T 

WRITE (6*5601 Till ,XCP,TANLLO(II 
DETER 1*TANLL0( I 1-TAND2 

XP1=( l./DETERII*(-TAND2* XCRdl +TANLLO! I)*(XP-YP* 

1TAND21) 

YP1=( 1. /DETER 1)*( -XCR( l )+XP-YP*T AND2) 

X4P=XP1 +SQRT ( ! X4PP-XP1) **2+(Y4PP-YPl)**2l*SIN(ATAN! (X4PP-XP l I /( Y4P 
IP— YP1 1 )— DELTA I 

X5P=X Pl+SORT! (X5PP-XP1 1**2+ (Y5PP-YP 11**21 *S INI ATAN( (X5PP-XP II /(Y5P 
1P-YP1 1 1 — DELTA ) 

B 1 = ( X4P— Sl*T ANL +CR/2. )/ (TANC-T ANL I 

B2 = ( X5P— SI*T ANO -CR/2. I / ( TANP-T ANO I 

B1=B1/S I 

B2=B2/SI 

83=1.00 

X4P=X4P/SI 

X5P=X5P/SI 

XP1=XP1/SI 

YP1=YP1/S I 

CR1=CR/SI 

B0T=1./SI 

TAPER1=TAPER 

CHDEXT=0. 

IF( B2.GE.0. I GO TO 2 

CHDEXT=CR l*( 1 .-T A PERI ) +.0001*! T ANP-TANO l+T AN0-TANL—B1*! TANC-T ANL I 
CR1=CR 1-CHDEXT 
TAPER 1=(X5P-X4P)/CR1 
B2=.000 1 

2 S=2.*(— 81**2* ( TANL— TANC ) /2. +82**2 *( TANO—TANP) /2.+B3*! CR1+B1*(TANL- 
1TANC) — B2* (TANO— TANP )+B3*(TAN0-T ANL I /2 . I +( 1.-B3 )/2.*(CRl+Bl* 

2( T ANL-T ANC )-B2*(T ANO-T ANP) +B3*(TAN0-TANL) ) I 
ARN=4./S 
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WRITE (6 ,553 ) 

WRITE! 6,554) X4P, X5P 

WRITE (6,563) ARN*CHI»LAP»P$I»0ME»B1, B2, TAPER1 »CR1 » BOT , CHOEXT 
WRI TE ( 6 ,566 ) 

WR ITE ( 6,565) XPI,YP1 
WRITE ( 6 ,568) 

4 CONTINUE 
GO TO l 

5 STOP 
END 
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WING IN LOW SWEEP POSITION 

ALL OF THESE DIMENSIONS ARE SCALED FOR A UNIT SEMISPAN IN THIS POSITION 
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WING IN LOW SWEEP POSITION 

ALL OF THESE DIMENSIONS ARE SCALED FOR A UNIT SEMISPAN IN THIS POSITION 
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